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Preface 


Situations  involving  both  intentional  and  inadvertent 
interference  in  communication  systems  are  very  common. 
Consequently,  there  has  been  substantial  interest  in  spread 
spectrum  modulation  as  a  method  of  mitigating  the  resulting 
impact  on  system  performance.  This  paper  investigates  three 
linear  spread  spectrum  formats  that  provide  the  inherent  ad¬ 
vantage  of  interference  rejection  with  the  unique  capability 
to  superimpose  desired  signals. 

I  would  like  to  express  my  appreciation  to  Captain 
Stanley  R.  Robinson,  my  thesis  advisor,  for  his  patient 
guidance  and  for  the  many  stimulating  discussions  offered 
throughout  the  course  of  this  work. 

This  project  would  not  have  been  possible  without  the 
help,  encouragement  and  understanding  of  my  wife,  Jackie. 

Her  unselfishness  for  the  past  two  years  has  permitted  me 
to  achieve  a  personal  goal  that  I  have  wanted  for  some  time. 

Frank  D.  Tiller 
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Abstract 

Spread  spectrum  techniques  can  produce  results  in  com¬ 
munication  and  navigation  systems  that  are  not  possible  with 
standard  signal  formats.  The  results  of  particular  interest 
are  the  ability  to  screen  messages  from  eavesdroppers,  high 
resolution  ranging  and  the  rejection  of  interference. 

Three  modulation  schemes  that  exploit  these  spread  spec¬ 
trum  properties  are  analyzed  in  this  thesis.  All  three  tech¬ 
niques  use  pulse-width  modulation  to  encode  sampled  analog 
voice  signals  onto  a  digital  phase  coded  carrier.  The  per¬ 
formance  analysis  assumes  transmission  in  the  presence  of 
white  Gaussian  noise  and  CW  jamming.  The  performance  analy¬ 
sis  uses  maximum  likelihood  estimation  to  define  the  sample 
parameter  estimate  and  associated  error  variance  when  prob¬ 
ability  density  functions  can  be  defined  for  the  receiver 
sample  outputs.  When  the  density  function  is  not  available 
a  second  moment  analysis  is  used  to  define  the  expected  sig- 
nal-to-noise  ratio. 

The  results  are  extended  to  include  the  introduction  of 
signal  dependent  code  noise  into  the  receiver  resulting 
from  the  degradation  of  the  orthogonal  code  word  structure. 
The  impact  of  code  noise  on  system  performance  is  seen  to 
grow  with  the  number  of  users  squared  and  the  received  sig¬ 
nal  energy. 
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PULSE-WIDTH  MODULATION 


OF 

SPREAD  SPECTRUM  CARRIERS 

I .  Introduction 

The  need  for  robust  communication  systems  has  caused  con¬ 
siderable  interest  within  the  last  several  years  in  spread 
spectrum  techniques.  The  application  of  spread  spectrum  mod¬ 
ulation  techniques  to  communication  systems  provides  capabil¬ 
ities  which  cannot  be  obtained  through  conventional  modulation 
formats.  The  primary  advantage  provided  by  spread  spectrum 
modulation  is  the  capability  to  reject  jamming  by  interfering 
signals  in  order  that  the  desired  signal  can  be  recovered. 

This  anti-jam  (AJ)  capability  is  of  great  significance  to  mil¬ 
itary  communications  as  it  provides  a  large  measure  of  protec¬ 
tion  from  active  countermeasures  by  an  enemy.  This  signal  dis¬ 
crimination  capability  also  allows  multiple  access  by  many  sig¬ 
nals  in  a  common  RF  channel. 

Spread  spectrum  modulation  results  in  a  wideband,  low-pow¬ 
er  density  signal  which  has  statistical  properties  somewhat  sim¬ 
ilar  to  random  noise.  As  a  result,  the  transmitted  signal  is 
not  readily  detected  or  recognized  by  conventional  surveillance 
techniques.  When  detected,  recovery  of  the  baseband  informa¬ 
tion  from  the  wideband  transmitted  signal  can  be  accomplished 
only  through  correlation  with  an  exact  replica  of  the 


transmitted  signal.  This  property  denies  the  casual  listener 
access  to  the  baseband  information.  Therefore,  spread  spectrum 
modulation  provides  message  privacy  and  selective  calling  cap¬ 
ability.  The  additional  inherent  navigation  capability  and 
high  multipath  tolerance  coupled  with  the  features  already  men¬ 
tioned  are  the  salient  reasons  for  employing  spread  spectrum 
modulation. 

Problem  Statement 

The  purpose  of  this  thesis  is  to  examine  and  evaluate 
three  direct  sequence  modulation  schemes  that  pulse-width 
modulate  the  pseudo-noise  codewords  (Dixon,  1976:  121)  with 
sampled  analog  voice  signals.  These  modulation  schemes  are 
antipodal  pulse-width  modulation,  bi-orthogonal  word  pulse- 
width  modulation  and  tri-orthogonal  word  pulse-width  modula¬ 
tion.  The  analysis  is  to  be  made  with  the  objective  of  pro¬ 
viding  a  mutual  system  performance  measure  for  the  three  en¬ 
coding  schemes  that  can  be  considered  when  choosing  a  modula¬ 
tion  format. 

Background 

There  are  four  basic  techniques  for  generating  spread 
spectrum  signals:  direct  sequence  modulation,  frequency  hop¬ 
ping,  pulsed-FM  or  chirp  and  time  hopping  (Dixon,  1976:  14). 
Direct  sequence  modulation,  as  the  name  implies,  is  generated 
when  a  binary  code  sequence  directly  modulates  the  radio  fre¬ 
quency  (RF)  carrier.  In  the  general  case  the  modulation  for¬ 
mat  may  be  AM,  FM  or  any  other  amplitude  or  angle  modulation 
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form.  The  modulation  techniques  examined  in  this  paper  all 
use  the  direct  sequence  or  pseudo-noise  format  to  bi-phase 
shift  key  (PSK)  the  RF  carrier. 

The  binary  code  sequences  to  be  studied  can  be  classified 
into  four  categories:  maximal  linear,  non-maximal  linear,  max¬ 
imal  nonlinear  and  nonmaximal  nonlinear  (Judge,  1962:81).  The 
maximal  linear  or  m  sequences  exhibit  a  two-valued  autocorre¬ 
lation  structure  for  whole  bit  increments  which  closely  approx¬ 
imates  that  of  a  random  binary  sequence.  The  maximal  nonlinear 
sequences  have  an  autocorrelation  which  is  a  function  of  the 
relative  time  shift  and  produce  a  many  valued  autocorrelation 
structure  similar  to  that  produced  by  the  linear  codes.  The 
nonmaximal  sequences  have  autocorrelation  properties  that  make 
them  inappropriate  for  most  practical  applications.  The  non¬ 
linear  maximal  sequences  are  used  in  most  practical  systems 
since  their  autocorrelation  structure  approximates  that  of  an 
m  sequence  and  they  offer  greatly  improved  privacy  due  to  the 
nonlinear  nature  of  the  code.  Linear  maximal  sequences  will 
be  assumed  in  this  paper  because  they  provide  an  upper  perform¬ 
ance  bound  for  non-linear  maximal  sequences  and  their  proper¬ 
ties  are  well  understood. 

Approach 

The  thesis  is  organized  into  six  main  sections.  The  first 
section  develops  the  notation  and  structure  used  throughout  the 
paper.  In  each  of  the  next  three  main  sections,  a  modulation 
>  format  is  analyzed  and  its  performance  is  predicted.  The 
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sequence  of  presentation  will  be  based  on  increasing  complex¬ 
ity  of  the  modulation  technique,  starting  with  antipodal  pulse- 
width  modulation  and  concluding  with  tri-orthogonal  word  pulse- 
width  modulation.  Performance  in  the  multi-user  environment 
and  code  noise  will  be  analyzed  in  the  sixth  section,  with  con¬ 
clusions  and  recommendations  presented  in  the  final  section. 

The  appendices  contain  a  summary  of  the  more  involved  issues 
necessary  for  the  analysis.  When  the  results  of  a  calculation 
or  derivation  are  included  in  the  appendix,  it  is  mentioned 
in  the  corresponding  section  of  the  text. 


II .  Signal  and  Channel  Model 


In  this  chapter,  the  basic  signal  and  channel  model  will 
be  defined  and  analyzed.  The  channel  model  will  apply  to  all 
pulse-width  modulation  formats  and  with  additional  definition 
the  general  signal  model  can  yield  any  of  the  three  modulation 
techniques. 

General  Signal  Model 

Maximal  linear  sequences  will  be  used  to  analyze  all  spread 
spectrum  modulation  formats  since  their  performance  provides  a 
very  tight  upper  bound  on  the  performance  of  more  practical 
nonlinear  maximal  sequences.  The  linear  sequences  also  approx¬ 
imate  very  closely  the  benchmark  characteristics  of  a  random 
binary  wave.  The  pulse-width  modulation  formats  require  at 
most  three  linear  pseudo-noise  sequences  that  are  mutually  orth¬ 
ogonal.  The  ith  code  waveform  will  be  written  as 

00 

Ci(t)  =  2  ajPT  (t_jTc)  (1) 

j—  C 

All  sequences  have  periods  T=NTc  where  N  is  the  number 

of  "chips"  per  codeword  and  T£  is  the  duration  of  a  chip  in 

seconds.  The  terms  a.  are  sequence  elements  taken  from 

J 

{+1,-1}  and  PT  (t)=l  for  0£t<T  and  zero  otherwise.  The 
Ac  c 

pseudonoise  sequences  are  thus  deterministic  waveforms  known 
only  to  the  transmitter  and  intended  receiver  that  repeat 
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every  T  seconds.  The  principle  characteristic  of  interest 
for  a  code  sequence  is  its  autocorrelation  function  which  is 
discussed  in  Appendix  A  and  the  Fourier  transform  of  the  auto¬ 
correlation  function,  its  power  spectral  density.  The  norm¬ 
alized  autocorrelation  function  of  a  linear  maximal  sequence 

—NT 

is  periodic  with  period  NT^  .  During  the  interval  c 


The  power  spectral  density  of  a  linear  maximal  sequence  is 
given  by 

S(f)  =  +  (ii2)  X)  sinoD‘(£-  w)  (3) 

k=-«  c 

kfO 

The  periodic  autocorrelation  function  and  discrete  power  spec¬ 
tral  density  of  a  pseudonoise  sequence  with  N=7  is  shown  in 
Figure  1  (Haykin,  1978:187). 

The  heart  of  the  general  transmitted  signal,  x(t)  ,  is  a 
generic  pseudonoise  waveform  c(t)  which  is  structured  to  sat¬ 
isfy  the  desired  pulse-width  modulation  format.  The  waveform 
structure  ranges  from  the  continuous  use  of  one  selected  wave¬ 
form  to  interleaving  segments  from  three  separate  codes.  The 
general  transmitted  signal  is  represented  by  a  binary  phase- 
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Fig  1.  Characteristics  of  Linear  Maximal  Sequence 

(a)  Autocorrelation  Function  (b)  Spectral  Density 

shift  keyed  (PSK)  format 

x(t)  =  /P  p(m,t)c(t)cos(coot)  (4) 

where 

wo  =  the  transmit  frequency  in  radian/sec.  (u)o=2irfo ) 

/ P  =  signal  amplitude 
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c(t)  =  the  pseudonoise  code  waveform  |c(t)|=l 
p(m,t)  =  pulse-width  modulation  |p(m,t)|=l 
mk  =  the  kth  sample  from  the  input  waveform  m(t) 
m  =  a  generic  sample  value 


This  transmitted  signal  is  limited  in  power  to  P/2  watts; 

therefore  the  total  energy  dissipated  in  T  seconds  cannot  ex- 
PT 

ceed  7^—  joules.  The  maximum  energy  associated  with  each  sam¬ 
ple  value  transmitted  is  given  by 


A  geometric  interpretation  of  the  transmitted  signal  as¬ 
sociated  with  each  modulation  scheme  analyzed  in  this  thesis  is 
available  through  signal  space  concepts  (Wozencraft,  1965:225). 
The  generic  pseudonoise  waveform,  c(t)  ,  is  composed  of  seg¬ 
ments  from  k  orthogonal  pseudonoise  waveforms  which  can  be 
thought  of  as  a  set  of  orthonormal  base  functions  for  the  sig¬ 
nal  space.  The  pulse-width  modulation,  p(m,t)  ,  is  used  to 
determine  the  coefficients  associated  with  each  bases  function 
in  the  k-dimensional  signal  space  and  therefore  determines  a 
unique  signal  vector  for  each  possible  sample  value.  The  sam¬ 
ple  values  are  obtained  from  sampling  the  normalized  in¬ 

put  waveform  m(t)  at  the  Nyquist  rate  of  ^  Hertz  as  shown 
in  Figure  2.  Because  the  sample  values  are  constrained  to 
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Normalized  Message  Sampling 


-IShkS+I  and  the  transmitter  is  power  limited,  the  trans¬ 
mitted  signal  will  lie  within  an  Euclidean  sphere  of  radius 
/PT/2  in  the  k  dimensional  signal  space.  The  receiver  will 
correlate  the  transmitted  signal  vector  with  the  bases  func¬ 
tions,  ,  to  obtain  the  transmitted  sample  value. 

In  an  actual  receiver  the  reconstructed  sample  value  can 
only  assume  N-l  discrete  values  corresponding  to  the  number 
of  chips  per  T  second  period.  When  the  input  waveform  m(t) 
is  assumed  to  be  distributed  uniformly,  the  mean  square  error 
due  to  quantization  is  given  by 


A/2 

\f 

-A/2 


x  dx  =  w  = 


12 (N-l) 5 


where 

A  =  width  of  one  quantization  level. 

N  =  number  of  chips  per  period. 

The  code  period  will  be  assumed  sufficiently  long  to  allow 
this  quantization  error  to  be  ignored  when  computing  the  per¬ 
formance  of  all  modulation  schemes. 


General  Channel  Model 

Using  a  tilde  to  denote  random  quantities,  transmission 
through  the  channel  results  in  the  additive  introduction  of 
the  noise  processes  Jj(t)  and  t  )cos(to0t-^)=^x(t  )cos(^0t )  + 

^y(t)sin(o)0t)  .  The  channel  noise  process  jj(t)  is  assumed 
to  be  a  zero  mean  white  Gaussian  process  with  two-sided 
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spectral  density  N0/2  .  The  second  noise  process  is  a  nar¬ 

row-band  jammer  with  random  uniform  phase  £  and  the  envel¬ 
ope  ^(t)  is  modeled  as  a  Rayleigh  process.  In  quadrature 


es  with  the  zero  means  and  the  same  autocorrelation  function 
given  by 

BJt(T)  =  J.sinc(Bt)  sl;llgT_l  (7) 


the  corresponding  power  spectral  density  is  given  by 

-  (r)n  (I)  <8> 

where 


H(t) 


I1*  1*1  -  * 

'o,  |t |  >  * 


O) 


The  channel  model  then  can  be  defined  as  a  discrete-time 
memoryless  additive  Gaussian  channel  with  an  average  power  con¬ 
straint.  The  channel  has  an  input  alphabet  and  channel 

output  alphabet  A  ,  both  equal  to  the  set  of  all  real  num- 

z 

bers  in  the  interval  [-1,+1]  .  If  x,,x2,x3,...are  the  inputs 
to  the  channel  at  times  1 , 2 , 3 , . . . then  the  corresponding  out¬ 
puts  z, ,z2 ,z3 , . . .are  given  by  zi=xx+yi  ’  where  yi»y2>y3>*** 
are  independent,  identically  distributed  normal  random  variables 
with  mean  zero  and  variance  a2  .  This  channel  is  often  depict¬ 
ed  as  in  Figure  3. 
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Y  :  Normal,  mean  zero 

a# 

variance  #* 


Fig  3.  The  Gaussian  Channel 


4.  Locus  of  Transmitted  Signal  Vector 


The  simplest  modulation  format  considered  will  use  one 

pseudonoise  code  or  bases  function  to  traverse  the  channel. 

The  signal  space  representation  of  this  modulation  scheme  is 

shown  in  Figure  4.  Allowing  mi  to  range  over  the  interval 

[~l,+l]  causes  the  tip  of  the  transmitted  vector  x  to 

—  m 

move  along  the  r,x  axis  from  -/E  to  +/E  .  The  modula¬ 
tion  technique  just  described  in  signal  space  is  antipodal 
pulse-width  modulation  and  will  be  analyzed  in  the  next  chap- 


Ill .  Antipodal  Pulse-Width  Modulation 


The  use  of  signal  formats  similar  to  antipodal  pulse- 
width  modulation  have  been  envisioned  for  more  than  ten  years 
as  an  efficient  method  of  communicating  between  terminals  in 
ranging  systems  using  pseudonoise  tracking  links  (Ward,  1967: 
69).  When  used  solely  to  provide  a  robust  communication  sys¬ 
tem,  antipodal  pulse-width  modulation  (PWM)  is  a  simple  direct 
sequence  method  of  transmitting  analog  speech  information. 

Modulation  Technique 

Antipodal  PWM  requires  only  one  code  waveform  Ci(t)  , 
assumed  to  be  a  pseudonoise  code  and  the  ability  to  employ  a 
phase  coherent  receiver.  The  transmitted  antipodal  PWM  sig¬ 
nal  x(t)  is  given  by 

x(t)  =  /F  d(t  )cos(u)0t )  (10) 

where  the  data  waveform,  d(t)  ,  is  defined  by 

d(t)  =  P(m,t)ci(t)  (11) 

and  the  pulse-width  modulation  is  represented  mathematically 
as 

P(m,t)  =  sgn'(t-kT-  |  -  (12) 

k=— oo 

with  sgn'(t)  defined  as  the  inverted  signum  function 


(-1,  t>0 
sgn'(t)  = I 

(+1,  t<0 

Combining  equations  (1)  and  (12)  results  in  the  expression 

00  00  /  y 

d(t>  "  I]  akPT  (t-1To>  S  sgn'(t-iT-  |  -  T-)  (13) 

k=-°°  c  i=-» 

i=fk/N] 

where  [k/N]  denotes  the  greatest  integer  value.  Physically 
each  analog  sample  is  mapped  into  a  code  period  by  the  pulse- 
width  modulation  signal  that  inverts  the  required  final  por¬ 
tion  of  the  period.  The  kth  analog  sample  value,  m^.  ,  de¬ 

termines  the  point  where  code  inversion  begins  during  the  in¬ 
terval  kT£t!(k+l)T  ,  which  corresponds  to  the  kth  code  per¬ 
iod.  The  start  of  both  the  code  period  and  PWM  signal  are 
synchronized  and  no  dead  time  is  allowed  between  PWM  signals. 
Figure  5  illustrates  typical  segments  from  the  code  and  PWM 
signals  before  they  are  mixed  and  applied  to  the  radio  fre¬ 
quency  (RF)  carrier  for  transmission.  The  general  correspon¬ 
dence  between  sample  values  and  the  alteration  produced  in 
the  transmitted  code  segment  is  shown  in  Table  I. 

Transmission  through  the  channel  described  in  Chapter  II 
results  in  a  received  signal  represented  by 

£(t)  =  /P  d(t  )cos(o)0t  )+Jj(t  )+jJ(t  )cos(u)0t-£)  (14) 

The  channel  disturbances  were  shown  to  be  approximately 
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TABLE  I 

Mapping  from  Sample  Value  to  Altered  Code  Segment 


m 

Upright  Interval 

Inverted  Interval 

-1.0 

0 

0  to  T 

-0.5 

0  to  T/4 

T/4  to  T 

0 

0  to  T/2 

T/2  to  T 

+0.5 

0  to  3/4T 

3/4T  to  T 

+  1.0 

0  to  T 

0 

Gaussian  and  represent  thermal  noise  and  a  narrow  band  jammer 
operating  at  the  receive  frequency. 

Receiver  Structure 

The  coherent  receiver's  structure  shown  in  Figure  6  rep¬ 
resents  a  charge  coupled  device  (CCD)  configured  as  an  analog 
correlator  (Melen,  19 77)-  The  sample  output  from  the  CCD  cor¬ 
relator  matched  to  ci(t)  can  be  written  as 


l  = 


2 


t )c \ (t ) cos (w  o t )dt 


(15) 


inserting  Jg(t)  into  the  integrand  gives 

T  T 

Tt  ~  2)/P  /  d(t  )c  i  (t  )cos2  ( to0 1  )dt  +  2  f  n  (t)Ci(t)cosU0t)dt+ 


(16) 


2  J*  ^(t )cos(to0t-£)ci (t )cos(to0t )dt 
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Fig.  6.  Coherent  Antipodal  Pulse-Width  Modulation  Receiver 
(Ziemer,  1976:321) 

The  narrow-band  noise  process  jj^(t)  which  appears  after  the 
RF  filter  in  the  receiver  can  be  represented  by  an  in-phase 
component  jjc(t)  and  a  quadrature  component  jjg(t)  •  Thus 

Jjf(t)  =  jjc(t  )cos(w0t )  -  jjs(t)sin(w0t)  (17) 

where  both  the  in-phase  component  jjc(t)  anc*  quadrature  com¬ 
ponent  jjg(t)  are  statistically  independent  Gaussian  random 
variables  with  zero-mean  and  variance  o2  (Papoulis,  1965: 
373). 

Using  the  trigonometric  identity  cos  ((jj0t  =  Hl+cos(2u)0t)) 
the  sample  output  can  be  written 
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T  T 

£  =  /P  /  d(t)cx(t)(l+cos(2u)0t))dt+  /  jjc(t)Cj  (t)(l+cos(2u)0t))dt- 

•'a  •'a 


2  J^Qs(t)ci(t)sin(<i)°t)  )  (  cos(aj0t|dt+2^  ci(t)cos(oj0t)  (18) 

[<^(t  )cos((i>  0t  )+,Jy(t  )sin((jj0t)jdt 


The  terms  with  [sin(oj0t )]  [cos(a)0t )]  in  the  integrand  re¬ 
duce  to  zero  since  the  sin  and  cos  are  approximately  orthogonal 
over  the  interval  0  to  T  Also  with  w0>>27r/T  the  double 
frequency  component  of  the  first  term  can  be  ignored.  The  as¬ 
sumption  is  valid  in  any  realistic  hardware  implementation  of 
the  correlation  receiver  since  the  frequency  response  is  such 
that  the  double  frequency  term  will  be  blocked.  The  correlator 
output  reduces  to 


d(t )cj (t )dt- 


,( t  )c  x ( t )dt+  /  Ci(t)J  (t)dt  (19) 


At  time  kT  the  sample  output  from  the  analog  correlator  can 
be  written  as 


2k  =  ^  mk  +  2  +  %  (20) 

where  m^  is  a  fixed  but  unknown  quantity  and  both  jj  and  ^ 

are  independent  Gaussian  random  variables.  The  density  func¬ 
tion  describing  the  random  variables  is  also  Gaussian 
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with  mean 

E{£k>  =  E(/F  T  mk>  +  E{jJ}  +  E{jj}  (21) 


'  /FT”k 

Since  the  random  variables  g  and  £  are  independent  with 
zero-mean,  the  variance  of  £k  is  found  from 


az  "  EV>  +  E^2> 


(22) 


■  /fv- 


3 )R  (a-8 )dad8+ 
c 


iFE°(‘- 

oo 


6) 


Rj(a-3)dad8 


Using  an  integral  transformation  where  x=ct-8  (Papoulis, 

1965:325)  and  the  approximation  that  R  (0-8)=  JU-  6( a-8) 

nc  z 

when  compared  to  the  code  auto-correlation  function,  the  var¬ 
iance  can  be  approximated  by 


0 


2 

z 


(23) 

(T-| x | )Rc(x)RJ(x)dr 


The  power  spectral  densities  of  the  code  waveform  and  nar¬ 
row-band  jammer  are  compared  in  Figure  7.  Examination  of  the 
spectral  densities  shows  B<<l/Tc  and  that  Sc(f)  is  approx¬ 
imately  flat  when  compared  to  the  jammers  power  spectral  den¬ 
sity.  Therefore,  R  (x)  can  be  written  as  an  impulse  function 
with  weight  T  .  Thus 

V 


20 


M+  y*  ( T-|t| )Tc6(T)J0sinc(Bx)dT  (24) 


The  sample  variance  is  then  approximately  given  by 


__ 2  „  N0T  J 0TJ 
az  '  + 


Signal  Processing  Performance 

The  output  of  the  analog  correlator  is  thus  a  Gaussian 


N0T  .  JoT2 


random  variable  with  mean  >/PTm  and  variance  2  N 

The  receiver  will  estimate  the  parameter  m  based  on  the  ob¬ 
served  sample  value  £  which  is  corrupted  by  additive  noise. 
Since  m  is  a  constant,  but  unknown  value,  a  logical  method 
(Van  Trees,  1968:63)  of  determining  m  is  to  find  inml  , 
the  value  of  m  that  maximizes  the  density  function  f(Z|m). 
The  quantity  m  ^  is  the  maximum  likelihood  estimate  or  the 
value  of  m  which  most  probably  gave  rise  to  the  observed 
quantity  Z  .  Thus  the  maximum  likelihood  estimate  of  the 
parameter  m  is  given  by 


*1  * 


The  receiver  performance  can  now  be  predicted  by  eval¬ 
uating  the  quality  of  the  estimator  structure  employed.  The 
first  measure  of  quality  is  the  expectation  of  the  estimate 


E{%1}  =  f  mmif(Z  |m)dZ  =  m 
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which  demonstrates  that  the  estimate  is  conditionally  unbiased. 
A  second  measure  of  quality  is  the  variance  of  the  estimation 
error 


Var{(“ml-m)  lm}=  E  {(%l“ra)2  lm} 


pfT 


(28) 


The  Cramdr-Rao  bound  (Van  Trees,  1968:66)  represents  a 
lower  bound  on  the  variance  of  any  unbiased  estimate  of  m 
and  is  given  by 


Var{ (m-m) | m}  > 


ainf(Zlm) 


9m 


11  -I 


- 1 


Solving  for  the  Cramer-Rao  bound  gives 


(29) 


- 

2 

Var{ (m-m) | m}  £ 

E 

Z-v'PTm 

a2 

z 

>  =  No/2  +  JqT/N 

PT2  PT 


(30) 


The  estimate  satisfies  the  bound  with  equality  and  is  termed 
an  efficient  estimate  of  m 

Since  a  closed  form  expression  for  the  density  function 
f(Z|m)  cannot  be  obtained  for  all  receiver  configurations,  a 
second  equivalent  measure  of  system  performance  will  be  de¬ 
fined  for  future  comparisons.  The  signal-to-noise  ratio  (S/N) 
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will  be  defined  as  (Wozencraft,  1965:588) 


S  _ 
N  " 


E2  {Z} 


(31) 


Thus  antipodal  PWM  results  in  a  signal-to-noise  ratio  given  by 


S  _  PT  m2 _ 

N  No/2  +  JoT/N 


Var{(mml-m) | m) 


(32) 


The  geometrical  ideas  associated  with  signal  space  give 
added  insight  into  the  performance  of  various  pulse-width  mod¬ 
ulation  formats.  As  discussed  in  the  last  chapter,  the  locus 
of  the  transmitted  vector  can  be  written  as 


X  .  m/E. 


(33) 


where  the  parameter  m  is  allowed  to  range  over  the  interval 
[l,+l]  .  The  transmitter  is  effectively  "stretching"  the  in¬ 

terval  [-1,+lJ  onto  the  larger  interval  [-/E,  +  v/E]  in  sig¬ 
nal  space.  The  "stretching"  is  uniform  in  the  sense  that 
(Wozencraft,  1965:611) 


dX 
— m 

dm 


=  SE 


for  all  m 


(34) 
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The  effect  of  maximum-likelihood  reception  is  to  undo  the 
"stretching"  performed  by  the  transmitter.  Therefore  the  sig¬ 
nal  "stretching"  is  related  to  the  receiver  performance.  De¬ 
fining  the  stretch  factor  S  by 


S 


(35) 


The  antipodal  PWM  signal-to-noise  ratio  can  now  be  written  as 


S  _  S2m2 
N  ~  N o  +2 J oT/N 


(36) 


where 


This  expression  for  S/N  shows  that  the  signal  energy  must  be 
increased  to  improve  performance  in  the  antipodal  PWM  system. 

The  next  logical  step  in  signal  space  would  be  a  signal 
format  with  two  bases  functions.  Figure  8  shows  the  signal 
space  representation  for  a  noncoherent  signal  format  with  two 
pseudonoise  codes  or  bases  functions.  The  locus  of  the  trans¬ 
mitted  signal  vector  can  be  written  as 

h, m  (tt)  i.  -  i.  <37> 

where  the  parameter  m  is  allowed  to  range  over  the  interval 
[-1,+1]  .  The  signal  format  just  described  in  signal  space 

is  bi-orthorgonal  word  PWM  and  will  be  analyzed  in  the  next 
chapter. 
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Fig  8.  Locus  of  Transmitted  Signal  Vector 
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V 


IV.  Bi-Orthogonal  Word  Pulse-Width  Modulation 

The  second  technique  for  information  transmission  in  a 
spread  spectrum  system  is  similar  to  antipodal  PWM.  The  ana¬ 
log  sample  value  controls  the  modulating  signal  as  before,  but 
now  varying  portions  of  two  orthogonal  codes  are  interleaved 
in  each  T  second  transmission  interval.  The  analysis  of 
both  a  coherent  and  noncoherent  system  will  be  given  in  this 
chapter. 


Modulation  Format 


The  generic  code  transmitted  in  bi-orthogonal  word  PWM, 
c(t)  is  composed  of  segments  from  two  orthogonal  pseudo¬ 
noise  waveforms  each^svith  T  second  periods.  Every  T  sec¬ 
ond  interval  of  the  generie^code  begins  with  Cj(t)  and  at 
the  point  determined  by  the  PWM  "signal ,  switches  to  c2(t) 
for  the  remaining  portion  of  the  intervqT^  The  duration  of 
each  sequence  portion  corresponds  to  the  analog  sample  itk  , 
such  that  when  nu  is  +1,  only  c^t)  is  sent^when  it  is 
-1,  only  c  2  ( t )  is  sent;  and  when  m.  is  0  ,  Ci(t)\is  sent 

i  i  \ 

\ 

for  the  first  half  of  the  interval  and  c2(t)  is  sent  for 
the  second.  There  is  no  dead  time  between  the  two  sequences. 
The  data  waveform  is  represented  mathematically  as 


d( 


t)  -  £  |c1(t)[u(t-kT)-u(t-kT-|  -  5^)] 


(38) 


m,  T 


+c2(t)  u( t-kT-  | - |-)-u(t-(k+l)T)]  | 
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where 


u(t)  = 


1,  t>0 
0,  t<0 


The  transmitted  spread  spectrum  signal  is  given  by 

x(t)  =  /P  d(t  )cos(co0t )  (39) 


Phase  Coherent  System 


After  propagating  through  the  channel,  the  resulting  re¬ 
ceived  signal  in  a  phase  coherent  bi-orthogonal  word  PWM  sys¬ 
tem  is  the  same  as  the  antipodal  PWM  system  with  an  alternate 


representation  for  d(t),  the  data  signal.  Thus  the  receiv¬ 
ed  signal  is  given  by 


Jg(t)  =  /P  d(t)cos(a)ot)+jj(t)+^(t)cos(u)ot-0)  (40) 

Receiver  Structure.  The  required  phase  coherent  receiver 

structure  is  shown  in  Figure  9.  The  output  from  the  ith  leg 

correlator  matched  to  codeword  c^(t)  can  be  written  as 

T 

|i  =  2  y*£( t)ci(t)cos(o)0t)dt  (41) 

o 


Now  working  with  ,  the  output  from  the  upper  leg  correla¬ 

tor,  matched  to  codeword  c^t)  and  expanding  the  integrand 


gives 


2/P  J  d(t)c1(t)cos2(w0t)dt  +2  J  nf  (t  )c  j  (t  )cos 
o  T  o 

( a)  0t  )dt+2  J  ?Jx(t)c1(t)cos2(a,0t)dt 
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2c2(  t)cosc4,t 


Fig  9.  Phase  Coherent  Bi-Orthogonal  Word  Pulse-Width 

Modulation  Receiver. 

Following  the  same  logic  presented  in  the  chapter  on  anti¬ 
podal  PWM  results  in  the  expression 

T  T  T 

£1  =  /P  /d(t)c1(t)dt+/^c(t)c1(t)dt  +  ftt  t)Cj(t)dt  (43) 

0  0  0 

Expanding  the  data  signal  d(t)  appearing  in  first  term  of 
equation  43  results  in 


« 


T 


b  - 

o 


c i (t )Ci (t )dt+/F 


c 

T 


T 

(t)c2(t)dt+  J  >jc(t)ci(t) 

o 


(44) 


dt  +  J c  i  (t  )^x(t  )dt 


The  second  term  in  equation  44  is  a  partial  cross-corre¬ 
lation  function  of  pseudonoise  codewords  Ci(t)  and  c2(t) 


over  the  final  segment  of  the  generic  code  period.  The  codes 
are  selected  to  be  approximately  orthogonal  when  the  crossn- 
correlation  occurs  over  a  complete  T  second  period,  however, 
cross-correlation  over  a  partial  code  period  will  not,  in  gen¬ 
eral,  be  zero.  The  partial  cross-correlation  function  then 
represents  an  additional  signal  dependent  interference  term 
for  the  receiver.  Chapter  VI  provides  an  analysis  of  this  sig¬ 
nal  dependent  noise  and  shows  that  the  partial  cross-correla¬ 
tion  in  channel  one  will  produce  an  additional  Gaussian  noise 
term  ^(m)  which  is  distributed  N(o,4'z(M)).  Then  at  time  KT, 
the  sample  output  from  the  correlator  in  the  upper  leg  can  be 
written  as 


h  -  +<M+£<m>  <45> 

where  m^  is  a  fixed  but  unknown  quantity,  ji  is  distributed 
N(o,NqT/2)  ,  £  is  distributed  N(o,J0T2/N)  and  ^(m)  is 

the  signal  dependent  noise. 

Following  the  same  logic  used  to  find  the  output  produced 
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in  channel  one,  the  output  from  the  correlator  matched  to 
c2(t)  in  the  lower  receiver  leg  is  given  by 

(t2)1  t 

^2  =  /P  c ! (t )c2 (t )dt+/F  j"  c2 (t )c2 (t )dt+ 

0  (l±m)T 


J  JJc(t)c2(t)dt+  f  c2(t)jJx(t)dt 


The  first  term  in  equation  46  is  the  partial  cross-correlation 

function  of  pseudonoise  codewords  Cj(t)  and  c2(t)  over  the 

initial  segment  of  the  generic  code  period.  Again  the  partial 

cross-correlation  function  represents  a  signal  dependent  noise 

to  the  receiver.  The  partial  cross-correlation  in  channel  two 

will  produce  an  additional  Gaussian  noise  term  £(m)  which  is 

distributed  N(o,A2(m))  .  The  noise  produced  in  each  receiver 

channel,  however,  has  a  symmetric  relationship  such  that  total 

cross-correlation  variance  a2  =  Y2 (m)+A2 (m)  is  not  signal  de 

c 

pendent.  Again  at  time  kT  the  sample  output  from  the  corre¬ 
lator  in  the  lower  leg  can  be  written  as 


{a  = 


+^+£(m) 


where  m^  is  the  fixed  but  unknown  parameter,  ri  is  distribut 
ed  N(o,N0T/2)  ,  is  distributed  N(o,J0T2/N)  and  £(m) 
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is  the  signal  dependent  noise  in  channel  two. 


The  final  ad  hoc  receiver  output  is  formed  by  subtracting 
the  output  of  receiver  channel  two  from  the  output  of  channel 
one  to  produce  the  Gaussian  random  variable  ^  . 

l  -  (48> 


The  expected  value  of  %  is  given  by 

E{Z}  =  E{Ij}-E{I2} 

=  v^Tm 


(49) 


Since  the  noises  in  each  channel's  sample  output  are  in¬ 
dependent  random  variables,  the  variance  of  Z  is  easily  com¬ 
puted  as 


a2_  = 


N0T  + 


2  J  0T2 
N 


+  4,z(m)  +  A2(m) 


=  N  nT  + 


2J0TZ 

N 


(50) 


+  o 


Signal  Processing  Performance.  The  observed  receiver 
sample  output  can  be  described  by  the  conditional  probability 
density  function  f(z|m)  as 

f(z|m)  =  — - —  exp 
/Sttcj 

z 

The  estimate  of  the  parameter  m  is  again  given  by 


( z-/FTm) : 


2‘z 


(51) 
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m  = 


z 

/PT 


(52) 


The  estimate  is  conditionally  unbiased  with  the  variance  of 
the  estimation  error  given  by 


Var{ (m-m) | m} 


PT2 


(53) 


NpT+2 J  0T2 /N+CTp 
PT2 


Neglecting  the  cross-correlation  noise,  the  variance  of 
the  estimation  error  is  double  that  seen  in  the  antipodal  PWM 
system.  Namely, 

Var{  (m-m)  |  m}  *  N°+2J°-T-Z.N.  (54) 

The  corresponding  signal-to-noise  ratio  reveals  the  same  3db 
degradation  in  system  performance 


S  PTm2 

N  N0+2J0T/N  v  ’ 

Phase  Noncoherent  System 

Noncoherent  bi-orthogonal  word  PWM  will  be  employed  when 
a  phase  offset  exists  between  the  received  signal  and  the  lo¬ 
cal  oscillator  in  the  receiver.  The  phase  offset  may  be  due 
to  doppler  shift,  oscillator  drift  and  other  common  distur¬ 
bances.  The  transmitted  signal  is  now  given  by 
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x(t)  =  /Pd(t  )cos(oj0t+(j)) 


(56) 


where  <f>  is  the  phase  difference  between  the  received  signal 
and  the  receiver  local  reference. 

Receiver  Structure.  The  noncoherent  bi-orthogonal  word 
PWM  receiver  is  shown  in  Figure  10.  The  complete  received  sig¬ 
nal  is  given  by 

Jt(t)  =  /Pd(t)cos(u)0t  +  <j>)+lj(t)+^(t)cos(<jOot-£)  (57) 

The  receiver  now  produces  inphase  and  quadrature  signal  com¬ 
ponents  at  the  output  of  the  correlator  matched  to  ck(t) 

The  correlator  outputs  are  given  by 
T 

=  2y'^(t)ck(t)cos(aj0t)dt,  k=lf  2  (58) 

o 


and 


Qk  =  2  /*£(t)ck(t)sin(u°t>dt>  k=l>2 


(59) 


Following  steps  similar  to  those  required  in  the  coher¬ 
ent  receiver  analysis  and  using  the  trignometric  identity 
sin2  ( co o t  )  =  z(  1-cos(2<a) 0t )  )  results  in  an  inphase  signal  given 
by 

T  T 

£k  =  /Pcosi}) j*  d(t)ck(t)dt+  J  nc(t)ck(t)dt+  (60) 

0  0 

T 

f ^x(t)ck(t)dt  ;  k-1,2 
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-Orthogonal  Word  Pulse-Width 
Receiver. 
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and  a  quadrature  signal  given  by 


£k  =  /Fsin$y’d(t)ck(t)dt-y')is(t)ck(t)dt+ 


(61) 


h 


(t)ck(t)dt  ;k=i,2 


Now  applying  the  definition  of  the  data  signal  d(t)  forces 
the  use  of  all  four  correlator  outputs.  The  equations  for 
outputs  produced  by  the  correlators  matched  to  codeword  one 

are 


b 


i/PCOS(J) 


(t)c1(t)dt 


0 


+ 


T 


/Pcos<(> 


cl(t)c2(t)dt 

T 


+ 


(62) 

T  T 

/ftc<t)Cl(t)dt  +  J C!(t)^x(t)dt 
o  0 


and 


/1+m^T 


IT"/'1'  T 

/Psin4>  f  c  i  (t  )c  i  ( t  )dt  +  ^Fsin<J(  /  ci(t)c2(t)dt 


(63: 


T  T 

f P„(t)Ci(t)dt  +  f J„(t)cx(t)dt 
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The  corresponding  sample  outputs  from  the  correlators  matched 
to  codeword  two  can  be  written  as 


C-?) 


£2  =  ^Tf-g — )cos<J>  +  +  4  +  cosij)j^(m) 


(68) 


and 


$2  =  /Pt(-k — ) sin<t>  -  JJ  +  %  +  sin<J>|(m) 


(69) 


where  as  in  the  phase  coherent  receiver  mk  is  a  fixed,  but 
unknown  quantity,  jj  is  distributed  N(o,N0T/2),  £  is  dis¬ 
tributed  N(o,J0T2/N)  ,  and  ^(m)  and  £(m)  are  the  signal 

dependent  cross-correlation  noises  distributed  N(o,'i'2(m))  and 
N(o,A2(m))  respectively. 

The  correlator  sample  outputs  are  therefore  independent 
Gaussian  random  variables  that  can  be  totally  characterized 
with  knowledge  of  their  mean  and  variance.  Table  II  provides 
the  mean  and  variance  associated  with  each  of  the  four  corre¬ 
lator  sample  outputs.  The  receiver  then  produces  random  var¬ 
iables  %  and  %  as  shown  in  Figure  10  which  are  defined  by 


S  =  /qx)2  +  (§i)2 


(70) 


%  - 


(71) 


The  probability  density  functions  for  %  and  %  cannot 
be  obtained  with  the  correlator  output  variances  defined  as 
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TABLE  II 


Correlator  Output  Statistics 


signal  dependent  and  unique  to  each  receiver  channel.  However 
Chapter  VI  shows  that  the  total  cross-correlation  noise 

(m)+A2  (m)  is  not  signal  dependent  and  in  this  case  the 
total  cross-correlation  noise  is 


a2  =  sin24>A2(m)+cos2<}>A2(m)+sin2<{>'t'2(m)+cos2<}!4'2(m) 
c 

=  y2(m)  +  A2 (m) 


(72) 


Then  to  obtain  a  closed  form  expression  for  the  conditional 
densities  f(x|m)  and  f(z|m)  ,  the  cross-correlation  noise 
in  each  channel  will  be  assumed  to  be  approximately  o2/4 
All  four  correlator  sample  outputs  now  have  an  approximate  var 
iance  given  by  N0T/2+J0T2 /N+o2/4 

The  random  variables  and  are  independent  and 
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jointly  normal,  with  the  means  given  in  Table  II  and  the  var¬ 
iance  in  each  channel  approximately 


(73) 


The  conditional  density  of  %  =  )2+(Qi )2is  given  by 

(Papoulis,  1965:498) 


f(x|m) 


(74) 


where 


and  I0(w)  is  the  modified  Bessel  function  of  order  zero  de¬ 
fined  in  Papoulis. 

The  conditional  density  of  g  =  /(£2 )2+($2 )i  is  obtain¬ 
ed  in  the  same  manner  and  is  given  by 


where 


f (z|m) 


2 


(75) 


Signal  Processing  Performance.  The  random  variables 
and  %  were  shown  to  have  a  Rice-Nakagami  or  Rician  condi¬ 
tional  probability  density  functions  in  the  last  section.  The 


final  ad  hoc  receiver  output  is  obtained  by  differencing  these 
two  Rician  random  variables  or  the  final  output  ^  is  given  by 


*  "  S  -  « 


(76) 


Given  the  ad  hoc  receiver  structure  shown  in  Figure  10,  the 
ideal  result  would  be  to  obtain  f(y|m)  ,  the  conditional 
probability  density  function  for  the  random  variable  y 
The  desired  density  could  be  obtained  from  the  convolution 
of  f(x|m)  and  f(z|m)  as 


fy(y|ffl) 


x(y-u)f 


z(-u)du 


(77) 


However,  this  integral  proved  formidable  and  this  route 
was  abandoned.  The  second  approach  was  to  use  the  character¬ 
istic  function  defined  as 


4>(o)) 


wxf (x)dx 


(78) 


41 


The  characteristic  function  for  the  random  variable  ^  is 
then  given  by  (Van  Trees,  1968:396) 


<i>Y(a>)  =  4>x(a))$z(-co) 


r  2  1 

exp 

jwn 

l-2jw6* 

(l^jooB2)* 


exp 


-jo)S2 

i+2jcoe2. 


(l+2ja>Ba)* 


(79) 


The  density  function  is  obtained  by  the  inversion  formula 


1 

r  “ 

j  exp 

r  *21 

jo on 

exp 

-jus2" 

.  1-2.] 

l+2j(jj32. 

2tt  J 

L  (l-2ju>e2)*  (l+2jioB2)i 

e-JU)ydu> 


(80) 


The  inversion  could  not  be  completed  and  the  conditional  den¬ 
sity  f(y|m)  was  not  obtained;  therefore  the  estimate  for 
m  and  its  associated  error  variance  could  not  be  found  for 
the  given  ad  hoc  noncoherent  receiver. 

The  alternate  signal-to-noise  ratio  analysis  requires 
the  first  two  moments  of  the  output  random  variable. 

The  nth  order  moments  are  found  for  the  random  variables  £ 
and  directly  to  be  (Middleton,  1960:414) 


(81) 


42 


where 


and  iFiCajBjix)  is  the  confluent  hypergeometric  function  dis 
cussed  in  Appendix  B. 

The  mean  is  therefore  given  by 


E(x} 


6iFi(-£;l;-a20) 


0e 


-n2/4B: 


(1+ 


(82) 


The  "large"  signal  approximation  for  the  expected  value  is  val 
id  when  n2>>82  is  given  by  (Papoulis,  1965:499) 

E{xjl.s.“  "(1+  fsr>  <83> 

Similarly  the  "small"  signal  approximation  is  given  by  (Pap¬ 
oulis,  1965:500) 

E(x>s.s.  “VFe<1+  (S4) 

Following  the  exact  same  logic  for  the  random  variable  ^  re¬ 
sults  in  a  "large"  signal  approximation  for  the  expected  val¬ 


ue  of 


and  a  "small"  signal  approximation  given  by 


E{z} 


s.s. 


(86) 


The  expected  value  of  the  random  variable  that  defines  the  re¬ 
ceiver  sample  value  output  is  found  from  E{y}  =  E{x}-E{z) 

The  "large"  signal  case  is  given  by 


-  S(lt  2iT> 


C'-s>+  T-  -  b 


(87) 


^PTm-  ^ —  (y-^-) 
/pT  1_m 


The  expected  value  of  the  random  variable  y  in  the  "small" 
signal  region  is  given  by 


(88) 


Comparing  the  results  obtained  here  for  the  noncoherent 
receiver  with  the  expected  value  of  r=v/PTm  obtained  for  the 
coherent  receiver  structure  indicates  some  distortion  effects 
are  introduced  by  the  noncoherent  receiver.  The  "large" 
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=  (0.313)Kr 

The  distortion  effects  are  shown  in  Figure  11  where  the 
expected  value  of  the  coherent  receiver  is  compared  to  the 
noncoherent  receiver  expected  value  approximations.  The  ap¬ 
proximations  are  given  as  a  function  of  K  ,  the  square 
root  of  the  signal  power  to  noise  power  ratio  and  show  the 
effect  of  variations  in  m 

The  second  central  moments  are  obtained  from  Middleton's 
result  (Equation  81)  and  reduce  to 
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E{x2 }  =  26z+n2 


(91) 


and 


E{z2 }  =  2g2+s2  (92) 

The  variance  is  the  required  parameter  for  the  signal-to-noise 
ratio  comparison.  The  variance  is  found  from  a2  =  E(x2)-E2(x) 
as 


o2=2g2+n: 

x 


(93) 


The  "large"  signal  approximation  for  the  mean  results  in  a 
corresponding  "large"  signal  approximation  for  the  variance 
given  by  (K=/PT/g) 


a 


2 

Xl  .  S  . 


*  6!|  r(lTi)!|  (94> 

-  3 2  1 1—  j^r(  1-m) 2  J 

where  the  last  equation  results  from  keeping  only  the  first 

°°  n  1 

two  terms  of  the  geometric  series  I  (-m)  =  — p - r-  and  im- 

n=0 

posing  the  condition  | m | < 1  on  the  modulation  sample  values. 

The  "small"  signal  approximation  for  the  variance  is  giv¬ 
en  by 
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-  26!+n!-  (§t)<4B!4t,!>2 


=  (2-  f)B2+d-  J)n2 


(95) 


*  $2((0.429)+(0. 054 )K2 ( 1+ra) 2 ) 


The  same  analysis  shows  that  the  random  variable  %  has  a 
"large"  signal  approximation  for  the  variance  given  by 


=  82(1-  ir( 1+m)2) 


(96) 


and  a  "small"  signal  approximation  for  the  variance  given  by 


a2  =  62 ( ( 0.429 )+(0. 054 )K2 ( 1-m) 2 ) 

zs.s. 


(97) 


The  total  variance  of  the  receiver  output  random  variable  v 


can  now  be  obtained.  Because  the  receiver  channels  produce 


independent  random  variables,  the  total  variance  is  given  by 


a2  =  a2  +  o2 


and  for  the  "large"  signal  region 


r  ~  202 ”  §T  (l+m2) 
1  •  s  • 


(99) 


and  the  "small"  signal  region  has  the  approximation 


s.s. 


*  (0.858)82+(0. 108) B2K2 ( 1+m2 ) 


(100) 
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The  coherent  receiver  variance  and  both  approximations  for  the 
noncoherent  receiver  variance  are  shown  in  Figure  12. 


! 


The  signal-to-noise  ratio  for  the  "small"  signal  region 
or  where  K  is  less  than  three  as  shown  in  Figures  11  and 
12  can  now  be  computed.  The  S/N  is  given  by 


E2{y} 


s.s. 


s.s. 


(101) 


_  1.94PT2  (^coherent 
8. 583 2 +1. 08PT2 ( 1+m2 ) 

,  PT2 

When  this  result  is  written  in  terms  of  K  =  ,  which  is 

the  signal  to  noise  ratio  without  the  modulation  parameter, 
the  impact  of  the  envelope  detectors  is  better  understood. 
The  signal-to-noise  is  then  expressed  by 


1.94  K2\N/coherent 

- z - 

8.58  +  1.08  K  (1+m2) 


1.94  K2 \ N  /coherent 
8.58  +  2.16  K2 


(102) 


Both  the  "small"  signal  and  "large"  signal  approximations  for 
the  signal-to-noise  ratio  are  plotted  in  Figure  13.  The  sig¬ 
nal-to-noise  ratio  plots  compared  to  the  coherent  signal-to- 
noise  ratio  isolated  the  distortion  produced  by  the  envel¬ 
oped  detectors. 
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Fig  12.  Noncoherent  Receiver  Variance  Approximations 
Compared  to  Coherent  Receiver  Variance 
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In  the  "large"  signal  region  or  where  K  is  greater  than 
three  the  signal-to-noise  ratio  is  given  by 


The  "large"  signal  region  signal-to-noise  ratio  converges  to 

the  same  (S/N)  obtained  in  the  coherent  receiver  structure. 

/PT 

In  the  region  of  operation  where  K=  —r—>5  ,  both  the  coher¬ 

es 

ent  and  noncoherent  bi-orthogonal  word  PWM  receivers  produce 
approximately  the  same  performance. 

The  fact  that  the  noncoherent  receiver's  S/N  approach¬ 
es  the  performance  obtained  with  the  coherent  receiver  can  be 
shown  by  an  alternate  method.  The  density  function  describ¬ 
ing  the  sample  outputs  from  either  envelope  detector  has  the 
form 

f(x|m)  =  f2exp[-  2^1]  Io(fr)  (104) 

with  ax>>82  this  density  function  can  be  approximated  by 
(Davenport,  1958:166) 
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(105) 


(x-g)2 

Tb5 


Hence,  when  the  magnitude  a  of  the  signal  component  is  large 
compared  to  3  (signal  standard  deviation)  and  when  £  is 
near  g  ,  the  probability  density  function  of  the  envelope  of 
the  sum  process  is  approximately  Gaussian.  Then  with  x  =  a 


f (x|m)=  — —  expf -  ]  (106) 

/2ttB  1  J 

The  envelope  detectors  produce  independent  random  var¬ 
iables  £  and  £  which  are  differenced  to  produce  the 
receiver  output  ^  can  now  be  approximated  (ax>>32)  by  a 
Gaussian  density 


f(y |m) 


- -  exp 

/2^/2W 

1  _  ( y-v'PTm) 2 

- -  exp  -  — - — 

y/^irv/2jT  l  2(232  ) 


y-(n-s)  I  2' 
2(232  ) 


(107) 


The  expected  value  and  variance  are  given  immediately  by 


E{y}1  =  /FTm 

X  •  s  • 


(108) 


and 


a 


2 


.S. 


=  232 


(109) 
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The  S/N  ratio  is  then  given  by 


S\  _  PT2m2 

N/ ,  B  ‘  2B2 


(110) 


This  S/N  result  is  the  same  obtained  by  both  the  original 
"large"  signal  approximation  and  the  coherent  bi-orthogonal 
word  receiver. 

The  signal-to-noise  ratio  can  again  be  defined  in  terms 

of  the  stretch  factor  introduced  in  the  last  chapter.  Since 

equal  increments  in  itk  correspond  to  equal  increments  in  the 

distance  measured  along  the  signal  locus  shown  in  Figure  8, 

the  stretch  factor  is  also  given  by  S=L/2  ,  where  L  is 

the  total  length  of  the  locus  traversed  by  x  as  m  in- 

—  m 

creases  from  -1  to  +1  .  The  length  in  the  bi-orthogonal 

PWM  format  is  /%E  ,  a  factor  of  /2  decrease  in  the  stretch 

factor  seen  in  the  antipodal  PWM  format.  The  signal-to-noise 
is  therefore 


S  =  (S2/2)m2 
N  No+2J0T/N 


(HI) 


The  signal  space  concept  of  signal  locus  length  then 
shows  the  same  3db  degradation  in  performance  and  points  to¬ 
ward  the  fact  that  the  signal  energy  must  be  increased  to  im¬ 
prove  the  bi-orthogonal  word  system.  A  logical  alternative 
would  seem  to  call  for  an  increase  in  the  dimension  of  the 
signal  space  to  produce  a  longer  signal  locus  and  corresponding 


performance  enhancement.  Figure  14  shows  a  three  dimensional 
signal  space  where  the  signal  locus  has  been  increased  to 
2  ✓SIT  which  projects  a  3db  improvement  in  system  performance 
over  the  antipodal  PWM  format.  Figure  14  is  the  signal  space 
representation  for  tri-orthogonal  word  PWM  that  will  be  dis¬ 
cussed  in  the  next  chapter. 
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Tri-Orthogonal  Word  Pulse-Width  Modulation 


The  final  technique  for  spread  spectrum  modulation  re¬ 
quires  segments  from  three  orthogonal  pseudonoise  waveforms. 
Tri-orthogonal  word  PWM  is  bi-orthogonal  word  PWM  with  one 
main  alteration  to  the  modulation  format.  The  new  feature 
seen  in  this  signal  format  is  the  selection  of  the  initial 
code  sequence  being  determined  by  the  polarity  of  itk  ,  the 
modulation  sample  value.  The  analysis  of  both  a  coherent  and 
noncoherent  system  will  be  given  in  this  chapter. 

Modulation  Format 

The  generic  code  transmitted  in  tri-orthogonal  word  PWM, 
c(t)  ,  is  composed  of  segments  from  three  orthogonal  pseudo¬ 
noise  waveforms  each  with  T  second  periods.  Every  T  sec¬ 
ond  interval  of  the  generic  code  begins  with  c3(t)  when  itk 
is  a  positive  quantity  or  begins  with  c3(t)  when  itk  is  a 
negative  quantity  and  at  the  point  determined  by  the  magnitude 
of  the  sample  value  switches  to  c2(t)  for  the  remaining  por¬ 
tion  of  the  interval.  The  duration  of  each  sequence  portion 
corresponds  to  the  analog  sample  rrK  such  that  when  iik  is 
+1,  only  Ci(t)  is  sent,  when  it  is  -1,  only  c3(t)  is 
sent,  and  when  m^  is  0  only  c2(t)  is  sent.  There  is  no 
dead  time  between  any  of  the  transmitted  sequence  segments. 

The  data  waveform  is  represented  mathematically  by 


57 


00 


d(t) 


]C  [o,<t) 

k=-oo 


u(t-kT)-u(t-kT 


+c2 (t ) 


u(t-kT- | 


m,. 


|  T)-u(t-(k+l)T) 


for  positive  mk  (0<mk-+l) 
and 


d(t) 


c3(t) 


u(t-kT)-u(t-kT 


+c2(t) 


u(t-kT-|mk|T)-u(t-(k+l)T 


(112) 


(113) 


for  negative  rak  (-l-mk<0) 

Figure  15  illustrates  how  positive  and  negative  sample 
values  control  the  interleaving  of  code  segments  to  form  the 
transmitted  generic  code  waveform.  The  transmitted  spread 
spectrum  signal  is  given 

x(t)  =  /P  d(t  )cos(oj0t )  (114) 

Phase  Coherent  System 

Using  the  same  channel  model  presented  in  Chapter  II  re¬ 
sults  in  a  received  signal  of  the  form 

J^t(t)  =  /P  d(t)cos(aj0t)+ji(t)+^(t)cos(u)ot-£)  (115) 

The  channel  has  introduced  both  a  white  Gaussian  noise  pro¬ 
cess  and  narrowband  jammer. 
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Fig  15.  Code  Sequence  Selection  for  Positive 
and  Negative  Sample  Values 


Receiver  Structure.  The  phase  coherent  receiver  struc¬ 
ture  for  tri-orthogonal  word  PWM  is  shown  in  Figure  16.  The 
output  from  the  ith  leg  correlator  matched  to  codeword 
ci(t)  can  be  written  as 

T 

=  2  j  Jg(t)ci(t)cos(u)0t)dt  (116) 

o 

During  each  T  second  transmission  interval,  either  the  out¬ 
put  ^  or  produces  just  a  noise  and  interference  com¬ 

ponent  with  no  signal  contribution.  Given  that  nu  is  posi¬ 
tive,  the  correlator  outputs  are  given  by 


i 


2c-j(  t)cosw0t 


irent  Tri-Orthogonal  Word  Pulse-Width 
lodulation  Receiver 
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*1  = /p  / 
0 

d(t)cj(t)dt  + 

J  J£c(t)ci(t)dt  + 

0 

J  jJx(t)ci(t)dt 
0 

T 

T 

T  (118) 

r» 

II 

* 

d(t)c2(t)dt  + 

f  jjc(t)c2(t)dt  + 

f  >Jx(t)c2(t)dt 

0 

Jo 

Jo 

and 


x  x 

h  =  y*c(t)c3(t)dt  +  Ax(t)c3(t)dt 

A 


(119) 


m|T 


^  J  ci(t)c3(t)dt  +  /P  y c2(t)c3(t)dt 


lm|T 


When  rtK  is  negative  there  is  a  symmetrical  result  where  the 
correlator  outputs  are  given  by 


l  ml  T 


&  =  y*Jec<t)ci(t)dt  +  A  (t)cl(t)dt  +  *  y c3(t)c1(t)dt 

A  A 


(120) 


+  </P  /  C2(t)C;(t)dt 

NT 

T  T 


^2  =  /P  d(t)c2(t)dt  +  jjc(t)c2(t)dt  +  ^x(t)c2(t)dt  (121) 


XX  X 

^3  =  /p y d(t)c3(t)dt  +  y ^c(t)c3(t)dt  +  y4(t)c3(t)dt  d22> 


Following  the  logic  used  to  evaluate  the  coherent  bi-orthogon- 
al  word  receiver  structure  produces  a  representation  for  pos¬ 
itive  and  negative  sample  values  that  has  a  symmetrical  struc¬ 
ture.  The  sample  outputs  at  time  kT  from  the  analog  corre¬ 
lators  can  be  written  as 
(l+m,  sgnm. 


/l+mksgnmk\ 

^ +  «  +  i  *  «<">) 

„  A_mksgn”lt\ 

*2  =  2 - 7Sgnmk  +  ft  +  4  + 


(123) 


(124) 


=  >1  +  4  +  4(m)  +  £(m) 


(125) 


where  is  a  positive  quantity. 


When  m^  is  a  negative  quantity  the  analog  correlators 


produce 


£1  “  $  +  4  +  + 


(126) 


/1-m  sgnmA 

£2  =  /pTf - - ^lsgnmk+  £  +  %  +  £(m) 

_  /1+mksgnm  \ 

Tf - *L- - —  )sgnmk+  jj  +  %  +  £(m) 


(127) 


(128) 


Where  in  both  cases  mk  is  a  fixed,  but  unknown  quantity,  jj 
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Signal  Processing  Performance.  The  symmetry  seen  in  the 
receiver  output  relations  insures  that  the  expected  value  and 
variance  of  the  random  variable  g  are  constant  over  time. 

The  observed  receiver  sample  output  can  therefore  be  describ¬ 
ed  by  the  conditional  probability  density  function  f(z|m)  as 


f ( Z | m) 


(Z-  /PTm) 2 
2A  2 


(132) 


The  maximum  likelihood  estimate  of  the  parameter 
by 


/\ 


m 


ml 


Z 

/PT 


m 


is  given 

(133) 


The  estimate  is  conditionally  unbiased  with  the  variance  of 
the  estimation  error  given  by 


Var{ (mml-m) | m} 


3NpT/2  +  3J0T2/N  + 


PT2 


(134) 


Neglecting  the  cross-correlation  noise  the  variance  of  the  es¬ 
timation  error  is  triple  that  seen  in  the  antipodal  PWM  sys¬ 
tem.  Namely 


Var  { (m  ^-m)  |  m} 


3N0/2  +  3JpT/N 
PT 


(135) 


The  corresponding  signal-to-noise  ratio  reveals  the  same 
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4.77db  degradation  in  system  performance. 


S  _  /1\  PTm2 
N  \  3  /  No_  J0T 
2  N 


(136) 


Phase  Noncoherent  System 

Noncoherent  tri-orthogonal  word  PWM  will  be  employed  when 
a  phase  offset  exists  between  the  received  signal  and  the  lo¬ 
cal  oscillator  in  the  receiver.  The  transmitted  signal  is 
then  given  by 

x(t)  =  /P  d(t  )cos(d)0t+4>)  (137) 

where  <J>  is  the  phase  difference  between  the  received  signal 
and  the  receiver  local  reference. 

Receiver  Structure.  The  noncoherent  tri-orthogonal 
word  PWM  receiver  is  shown  in  Figure  17.  The  complete  re¬ 
ceived  signal  is  given  by 

JSJ(t)  =  /Pd(t  )cos(uj  0t+4>  )+n(t  )+jJ(t  )cos(u)0t-J)  (138) 

The  receiver  now  produces  inphase  and  quadrature  signal  com¬ 
ponents  at  the  output  of  the  correlator  matched  to  c^t) 

The  correlator  outputs  are  given  by 

T 

Jk  -  2 /  8(t)ck(t)cos(u0t)dt,  k=l,2,3  (139) 

*'o 

and 

T 

$k  =  2/s(t)c k(t)sin(w0t)dt,  k=l>2,3  (140) 
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Filter 


Fig  17.  Phase  Noncoherent  Tri-Orthogonal  Word  Pulse-Width 

Modulation  Receiver 
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Following  the  steps  outlined  in  the  noncoherent  bi-ortho- 
gonal  word  receiver  analysis  results  in  an  inphase  signal 
given  by 


=  /Pcos*})  /  d(t)ck(t)dt 

**  ft 


1 

*/ 


$c(t)ck(t)dt 


1 


5x(t)ck(t)dt 


;  k-1,2,3 


(141) 


and  a  quadrature  signal  given  by 


T  T  T 

§k  =  /Fsin<f)  /  d(t)ck(t)dt-  /  jjs(t  ,ck(t  )dt+  J  ^y(t)ck(t)dt 

**  *  'n  *  ft 


(142) 


;  k-1,2,3 

The  representation  for  positive  and  negative  sample  val¬ 
ues  has  the  same  symmetrical  character  seen  in  the  analysis 
of  the  coherent  receiver  structure.  When  mk  is  positive 
the  sample  value  outputs  at  time  kT  of  the  correlators  mat¬ 
ched  to  Ci(t)  are  given  by 


/l+mksgnmk\ 

=  /PTcos<{>( - j  sgnmk+jj+,J+cos<j>^(m) 


and 


/l+m  sgnm.\ 

Qi  =  ^Tsin^l - - - ■)  sgnmk-^+J+sin(j)^(m) 


(143) 


(144) 
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1 


Ia 


The  sample  outputs  from  the  correlators  matched  to  codeword 
two  can  be  written  as 


_  /1-m  sgnnuX 

£2  =  /PTcos<{>/ - I  sgnm^+^+^+coscfiCC 


m) 


(145) 


and 

/1-m.sgnm. 

j^2  =  /PTsin<j>( - J  sgnmj^-jj+^+sin^CCm) 


■) 


(146) 


The  corresponding  sample  outputs  with  positive  mk  for  chan¬ 


nel  three  are  given  by 


£3  =  £  +  +  cos^Cm)  +  cos^(m) 


(147) 


and 


§3  =  -JJ  +  £  +  sin<t>£(m)  +  sin<K(m) 


(148) 


When  m^  is  a  negative  quantity  the  outputs  from  channel  one 
and  channel  three  are  interchanged.  The  sample  outputs  from 
the  correlators  matched  to  codeword  one  are  thus 


li  =  %  +  fl  +  cos^(m)  +  cos<K(m) 


and 


$1  =  -JJ  +  +  sin<f>£(m)  +  sin<J)f  (m) 


(149) 


(150) 


The  correlator  matched  to  codeword  two  produce  outputs  defined 
by 
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_  /l-mksgnmk\ 

^2  =  /PTcos4>( - !  sgnmk  +  jj  +  ^  +  cos<f>£;(m) 


(151) 


and 


_  /1-mkSgnmk\ 

=  /PT  sin<J,( - - ±J 


sgnmk  -  jj  +  i[  +  sin<|>£  (m)  (152) 


The  corresponding  outputs  from  the  correlators  in  channel 
three  are  given  by 


/l+m  sgnmA 

^3  =  /PT  cos<f>l - )sgnmk  +  ji  +  +  cos^(m) 

and 

_  /1+m.sgnm  \ 

$3  =/PT  sin<{>f - - - J  sgnmk  -  jj  +  %  +  sin<t>£(m) 


(153) 


(154) 


Where  in  both  cases  mk  is  a  fixed  but  unknown  quantity,  jj 
is  distributed  N(o,N0T/2)  ,  jJ  is  distributed  N(o,*^T  ■)  , 

£(m)  is  distributed  N(o,H'2(m))  ,  and  £(m)  is  the  comple¬ 

mentary  signal  dependent  noise  distributed  N(o,A2(m)). 

The  correlator  sample  outputs  are  thus  Gaussian  random 
variables  that  can  be  totally  characterized  with  knowledge  of 
their  mean  and  variance.  Table  III  provides  the  mean  and  var¬ 
iance  associated  with  each  correlator  output  and  sample  para¬ 
meter  polarity.  The  receiver  then  uses  envelope  detectors  to 
produce  intermediate  outputs  ^  ^  and  %  as  shown  in 

Figure  17.  The  envelope  detector  outputs  are  defined  by 

3  =  +  CQVF  (155) 
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Table  III 


Correlator  Output  Statistics 


Correlator  Output 

Mean 

Variance 

II 

/PT 

1+m] 

L2  J 

COS(J) 

N0T 

2 

+ 

J0T2 

N 

+cos2  4>4'2  (m) 

Qi 

v'PT 

1+m  | 

-  2  J 

sin<}> 

N0T 

2 

+ 

J0T2 

N 

+sin2<j>4'2  (m) 

s 

w 

> 

I2 

y'PT 

1-m] 

-  2  J 

COS(() 

N0T 

2 

+ 

J0T2 

N 

+cos2  (()A2  (m) 

HH 

w 

o 

Qz 

/PT 

”l-m* 

2 

sin<}) 

N0T 

2 

+ 

J0T2 

N 

+sin24>A2  (m) 

ft 

I3 

0 

N0T 

2 

+ 

J0T2 

N 

+cos24>a2 

c 

Q3 

0 

N0T 

2 

+ 

J0T2 

N 

+sin  z<boz 
c 

Ii 

0 

N0T 

2 

+ 

J0T2 

N 

+cos2d>o2 

c 

Qi 

0 

N0T 

2 

+ 

J0T2 

N 

+sin2<j>a2 

w 

> 

i2 

/PT 

1-m 

L2  J 

COS((l 

N0T 

2 

+ 

J0T2 

N 

+cos2  c}> A 2  (m) 

Eh 

< 

O 

w 

53 

Q2 

/PT 

1-m 

2 

sin4> 

N0T 

2 

+ 

J0T2 

N 

+sin2  <J>A2  (m) 

l3 

/PT 

1+m 

.  2  J 

cos<(> 

N0T 

2 

+ 

J0T2 

N 

+cos2  (m) 

q3 

/PT 

fl+m 

L2  . 

sin<{> 

N0T 

2 

+ 

J0T2 

N 

+sin2<J>4'2(m) 

y 


(156) 


(157) 


The  probability  density  functions  for  £  ,  ^  and  z  cannot 
be  obtained  with  the  correlator  output  variances  defined  as 
signal  dependent  and  unique  to  each  receiver  channel.  Chap¬ 
ter  VI  shows  that  the  total  cross-correlation  noise  a2=A2 Cm)  +  4'2 

c 

(m)  is  not  signal  dependent  and  as  shown  in  Chapter  IV  each 
envelope  detector  channel  produces  approximately  a2/4  cross¬ 
correlation  noise.  Using  this  approximation,  the  correlator 

outputs  with  non-zero  means  all  have  a  variance  equal  to 
9  a2 

NoT  JoT  c  ,  ,  ..  . 

— 7j—  +  — and  the  correlator  outputs  with  zero  means 

NoT  <J  T2 

have  a  variance  equal  to  — ^ .  The  probability 
density  functions  can  now  be  defined  using  the  approximate 
channel  variances. 

Signal  Processing  Performance.  The  final  receiver  out¬ 
put  y/  is  formed  from  a  linear  combination  of  the  envelope 
detector  outputs.  The  receiver  output  is  defined  by 

«  =  S  +  £  "  K  (158) 

The  symmetry  seen  in  the  envelope  detector  outputs  insures 
that  the  expected  value  and  variance  of  the  random  variable 
yt  are  constant  over  time  since  interchanging  the  roles  of 
^  and  %  has  no  impact  on  jv.  Therefore,  throughout  the 
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remainder  of  this  section  ^  will  be  assumed  to  be  the  envel¬ 
ope  detector  output  with  zero  mean. 

The  random  variables  and  ^  are  independent  and 

jointly  normal  with  means  given  in  Table  III  and  the  variance 
in  each  channel  approximately 


2  _  N0T  J0T2  °c 

2  N  4 


(159) 


The  conditional  density  of  x=/(^i ) 2  +  ( Q  i ) 2  is  given  by 
(Papoulis,  1965:498) 

f(x|m)  =  jr  exp  X- +^-  -  j  I0fp-^U(x)  (160) 
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The  conditional  density  g  =  /(£3 ) 2  +  (Q3 ) 2  is  obtained  when 
both  and  j^3  are  zero  mean  as  a  simplification  of  the 

usual  Rician  density  function  (Raemer,  1969:85).  The  result¬ 
ing  envelope  is  distributed  Rayleigh  with  the  probability  den¬ 
sity  function  given  by 


f(z|m) 


exp 


(162) 


2  (  6  2  + 


U(Z) 


Thus,  the  random  variables  £  and  ^  have  a  Rician  condi¬ 
tional  density  function  while  the  third  channel  ^  has  a 
Rayleigh  density  function.  The  final  receiver  output  is 
formed  by 


%  =  S  +  ?  “  X,  (163> 

The  density  function  f(wim)  could  not  be  obtained  as 
shown  in  Chapter  IV.  Thus,  an  estimate  for  m  and  its  as¬ 
sociated  error  variance  could  not  be  found  for  the  given  ad 
hoc  noncoherent  receiver.  The  alternate  signal-to-noise  ra¬ 
tio  requires  only  the  first  two  moments  of  the  output  random 
variable.  The  expected  value  is  found  from 

E{w}  =  E{x)  +  E{z)  .-  E{y }  (164) 

The  results  derived  in  Chapter  IV  and  E{Z}  = ^ ^( 82+c2 /2 ) 
immediately  provide  the  expected  value  of  the  receive  ’  out¬ 
put.  Thus,  in  the  "large"  signal  region  (PT2m2>>82) 
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and  in  the  corresponding  "small”  signal  region 


Etw}s.s.  If  ".+^5  [B*  +  %-]  (166) 


Comparing  the  results  obtained  here  for  the  noncoherent  re¬ 
ceiver  with  an  expected  value  of  r=/PTm  obtained  for  the 
coherent  receiver  structure  indicates  some  distortion  effects 
are  introduced  by  the  noncoherent  receiver.  Neglecting  the 

cross-correlation  noise,  the  "large"  signal  approximation 

/PT 

can  be  written  with  K=  — as 

p 


E{w}l.s.  “  r“  W  ( 1+m2  )+  [5] 


(167) 


The  "small"  signal  approximation  can  also  be  written  in  terms 
of  the  coherent  receivers  expected  value 


Etw,s.s.  “  °-313Kr  +  k  VI  (168> 

The  distortion  effects  are  shown  in  Figure  18,  where  the  ex¬ 
pected  value  of  the  coherent  receiver  is  compared  to  the  non¬ 
coherent  receiver  expected  value  approximations.  The  approx¬ 
imations  are  plotted  as  a  function  of  K  ,  the  square  root  of 
the  signal  power  to  noise  power  ratio.  The  distortion  caused 
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I 


Fig  18.  Noncoherent  Receiver  Expected  Value  Approximations 
Compared  to  Coherent  Receiver  Expected  Value 

\  ^  O  / 


6 
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by  the  direct  feed-through  of  noise  in  one  channel  is  very 
severe  when  the  signal-to-noise  ratio  is  small.  The  receiv¬ 
er  should  perform  poorly  in  a  high  noise  environment. 

The  corresponding  "large"  signal  approximation  for  the 
receiver  output  variance  is  given  by 

oj  =  2g2 -  ( 1+m2 )  +  (2-  £)B2  (169) 

1  •  s . 


and  the  "small"  signal  region  has  the  approximation 

(170) 

o2  =  ( 0 . 858 ) $z  +  (0. 108) 6 2K2 (1+m2 )+(2-  £)g2 

s.s.  * 


The  coherent  receiver  variance  and  both  approximations  for 
the  noncoherent  receiver  variance  are  shown  in  Figure  19. 

The  signal-to-noise  ratio  for  the  "small"  signal  region 
or  where  K  is  less  than  three  as  shown  in  Figures  18  and  19 
can  now  be  computed.  The  S/N  is  given  by 


E2{w} 


s.s. 


w  s.s. 


57T  |  TT7T  m  +  ej 


12. 9B*  +  1. 08f52K2(l+mz) 


(171) 


=  15n( /PT/4g+g//PTm)2 
12 . 9+1 . 08PT2 ( 1+m2 )/g2 


(!) 


coherent 


In  the  "large"  signal  region  or  where  K  is  greater  than  three 

V 
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the  signal-to-noise  ratio  is  given  by 


/PTm-(B2//PT)(l+m2)-V2  8 

2 .43B2-( $4 /PT2 ) ( 1+m2 ) 


3 

1_  FFm  (1+m  )+^" 

B  ' 
i/FTm. 

2 

2.43  -  ^  (1+m2)] 

\N/ 

coherent 


(172) 


The  "large"  signal  region  signal-to-noise  ratio  does  not  con¬ 
verge  to  the  S/N  obtained  for  the  coherent  structure  due  to 
the  distortion  introduced  by  the  third  channel  feed  through. 

Both  the  "large"  and  "small"  signal  approximations  for  the  sig¬ 
nal-to-noise  ratio  are  plotted  in  Figure  20.  The  signal-to- 
noise  ratio  approximation  plots  compared  to  the  coherent  signal- 
to-noise  ratio  isolates  the  distortion  produced  by  the  envel¬ 
ope  detectors.  The  third  curve  shows  the  true  performance  when 
the  distortion  terms  are  removed. 

The  signal  space  model  also  predicts  the  same  4.77db  de¬ 
gradation  in  system  performance  seen  in  the  coherent  tri-orth- 
ogonal  word  receiver  when  viewed  correctly.  The  stretch  fac¬ 
tor  for  any  sample  value  is  /E/2  rather  than  the  expected 
/2E  value.  Also  an  additional  noise  term  is  introduced 
through  the  third  bases  function  to  give 


■S  =  _ s! _ 

N  V3/No+2J0T/N 


(173) 
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(ti/s)  % 


Fig  20.  Noncoherent  (S/N)  Approximations  Versus 
Coherent  (S/N) 
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Thus  the  additional  noise  seen  in  the  tri-orthogonal  word  for 
mat  degrades  the  expected  system  performance  below  that  pro¬ 
vided  by  the  earlier  modulation  formats.  Because  of  the  sev¬ 
ere  impact  of  additional  noise  sources  in  any  modulation  for¬ 
mat,  the  next  chapter  will  examine  signal  dependent  noise  pro 
duced  in  the  multiple  code  word  receivers. 


VI.  Code  Noise 


All  multiple  codeword  pulse-width  modulation  schemes  pro¬ 
duce  a  partial  cross-correlation  function  in  addition  to  the 
desired  signal.  The  code  cross-correlation  function  will  rep¬ 
resent  a  signal  dependent  noise  source  to  the  receiver.  This 
chapter  will  examine  code  cross-correlation  noise  in  both  the 
single  and  multiple  user  environment.  The  random  binary  wave 
discussed  in  Appendix  A  will  be  used  to  set  bounds  on  the  per¬ 
formance  of  any  deterministic  code  sequence  used  in  an  actual 
communication  system. 


Self  Noise 


The  coherent  bi -orthogonal  word  PWM  system  will  be  used 
to  develope  the  partial  cross-correlation  noise  results  that 
apply  to  all  multiple  codeword  modulation  formats.  The  out¬ 
put  of  just  the  signal  component  from  the  analog  correlators 
was  given  in  Chapter  IV  as 

jl+m|T  „T  (174) 


Ii  =  /P 


\72-f  r 

I  ci(t)cj(t)dt  +  /P  I  ci 

3 


(t  )c2  (t  )dt 


and 


/P 


C 

, (t )c2 (t  )dt  +  I  c2(t)c2(t 


(175) 

)dt 
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Both  receiver  channels  produce  a  desired  signal  term  and  an 
unwanted  signal  dependent  noise  term.  Figure  21  illustrates 
for  a  typical  transmitted  code  sequence  the  formation  of  both 
the  desired  term  and  cross-correlation  noise  in  the  receiver 
channels.  The  correlator  outputs  reduce  to 


♦  ?(m) 


(176) 


and 


(177) 


The  output  from  channel  I  is  then  a  signal  term  and  par¬ 
tial  cross-correlation  function.  The  partial  cross-correla¬ 
tion  function  is  given  by 

fT 

^(m)  =  /P  I  c2(t)ci(t)dt  (178) 


Assuming  that  both  Cj(t)  and  c2(t)  are  random  binary  waves 
composed  of  equally  likely  pulses  of  amplitude  ±1  and  dur¬ 
ation  Tc  seconds,  their  temporally  aligned  product  Cj(t) 
c2(t)  will  produce  another  random  binary  wave  Cs(t)  with  the 
same  characteristics. 

The  partial  cross-correlation  function  is  then  expressed  by 

r 

£(m)  =  /P  /  c  3 ( t )dt  (179) 

{¥}T 
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Fig  21.  Origin  of  Signal  Dependent  Noise  in  Multiple 

Codeword  Systems 


The  first  two  moments  of  the  signal  dependent  noise  can  now  be 
computed.  The  expected  value  of  ^(m)  is  given  by 


E{e(m) } 


E{c3 (t ) }dt 


(180) 


=  0 
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The  variance  of  the  signal  dependent  noise  is  then  found  from 

T  T 

E{e2 (m) }  =  P  j  J  E|c3(ti  )c3(t2  )Jdtidt2 
{^}T  {^}T 


(181) 


/  / 

i"i  (j?). 


Rc(t1-t2)dtldt2 


Defining  the  variable  x~t1-t2  ,  the  double  integral  can  be 

transformed  into  a  single  integral  in  t  as  shown  in  Fig¬ 
ure  22  (Papoulis,  1965:335).  The  variance  is  then 


A 

E{e2 (m) }  =  P  J  (A-|x| )Rc(T)dx 


(182) 


where  R  (x)  is  the  autocorrelation  function  of  the  random 
c 

1 —in 

binary  wave  defined  in  Appendix  A  and  A=T(— ^— )  . 

Substituting  the  expression  for  RC(T)  into  the  inte¬ 


grand  results  in 


E{e2 (m) }  = 


/ 


1  -  |  x  1  ][l-  -^JdT  (183) 


i 


i 


Fig  22.  Double  Integral  Transformation 


■ 2P/ 


The  second  channel  also  produces  a  partial  cross-correla¬ 
tion  term  defined  in  Equation  175  and  given  by 
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c3(t)dt 


(184) 


£(m)  =  /P 


/ 


The  expected  value  is  found  to  be 

{ttJ* 


E{£(m)}  =  /P  |  E(c3(t)}dt 


(185) 


The  variance  of  the  signal  dependent  noise  in  the  second  chan¬ 
nel  is  found  from 

EU2(m)}  =  P  /  /  E{c3(t1)c3(t2|dt1dt2 


{t=}T  {¥}t 

•/  / 


(186) 


R  (t i — 1 2 )dt idt2 


Using  the  same  integral  transformation  required  for  channel  I 
results  in  a  single  integral  expression.  The  variance  in 


channel  II  is  given  by 


E{£2  (m) }  =  P  f  [T|ij2|-|T|  J  rc(t)  dx 
-{^}T 
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PTT 
_ c 

2 


PT  2 

( 1+m)- 


Defining  the  variance  in  channel  I  as  f2(m)  and  the  variance 
in  channel  II  as  A2(m)  allows  the  total  cross-correlation 
variance  o2  to  be  defined  since  the  channel  outputs  are  tem- 

C 

porally  juxtaposed  and  therefore  are  independent  random  var¬ 
iables.  The  total  cross-correlation  variance  is  not  signal  de¬ 
pendent  and  is  defined  by 


a2  =  A2(m)  +  ^(m) 

c 


=  PTT 

c 


2PT2 

c 


3 


(188) 


The  partial  cross-correlation  noise  is  formed  by  adding  (in¬ 
tegrating)  the  contributions  of  a  large  number  of  random 
"chips"  in  the  random  binary  waveform.  The  central-limit 
theorem  then  insures  that  the  resulting  cross-correlation 
noise  will  be  approximately  Gaussian. 

Multiple  User  Code  Noise 

The  results  derived  in  the  last  section  can  be  extended 
to  multiple  users  operating  in  a  conferencing  mode.  The 
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correlator  outputs  in  the  coherent  bi-orthogonal  word  receiv¬ 
er  now  contain  the  superposition  of  the  K  users  in  the  confer¬ 
ence.  The  correlator  outputs  are  given  by 


and 


- 


XX  j 

i=l  * 

K  { 
I>J 


(189) 


(190) 


The  expected  value  of  the  partial  cross-correlation  noise  is 
given  by 


K 

y]  E{ei(m)} 
i=l 


=  0 


and 


K 

^  E{Ci(m)} 
i=l 


(191) 


=  0 


(192) 


The  partial  cross-correlation  variance  produced  in 


channel  I  of  the  receiver  is  found  from 


.  K  K  .  K  K 

{2  2<i<n,)*N('n)}  ■  ££ 

'  i=l  N=1  '  i=l  N=1 


E{ei(m)eN(m)} 


(193) 


K  K  T  T 

=  J  f 

i=l  N=1 


R  (t  j-t2 )dt idt2 


The  sample  values  mj ,m2 ,m3 . . .m. ,m  .  ..m„  can  be  assumed  to 
occur  as  increasing  algebraic  numbers  in  the  interval  [-l,+l] 
without  loss  of  generality.  This  result  occurs  since  any 
other  ordering  of  the  K  sample  values  is  only  a  permutation 
of  this  ordering  and  all  permutations  sum  to  the  same 
quantity.  Now  assuming  that  P^P^P  and  m^n^  ,  the  dou¬ 
ble  integral  given  in  equation  193  can  be  transformed  to 
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■■  I'll  ^ 1 

- 1*,  <2n-1>  ~irL  <1-^)  -  -r- 


The  partial  cross-correlation  variance  produced  in  chan¬ 
nel  II  of  the  receiver  is  found  using  the  same  approximations. 


The  variance  is  then  given  by 


K  K 


,K  k  .  K  K  Tji^jTja-a( 

-EE¥.  f  ( v 

*i=l  n=l  '  i=l  n=l  J  J 


(195) 


R  (t i-t2 )dt  jdt2 


The  sample  values  again  occur  as  increasing  algebraic 
numbers  in  the  interval  [-l,+l]  where  m^<m^  and  P^~P^=P  . 
The  double  integral  in  equation  195  can  now  be  transformed  to 

Efe  S5i(mHn(m)}  ‘  Z)(2K-2i+1)P  /  [T<ii^HTi][l-4L]d* 

li=l  n=l  '  i=l  »/_  i  C 


(196) 


K  p  2 

=  ^(2K-2i+l)^nr£(l+mi)- 


The  results  for  the  multiple  user  code  noise  reduce  to  those 
obtained  for  the  self  noise  produced  by  a  single  user.  The 
total  cross-cc  ^relation  noise  produced  in  the  multiple  user 
system  is  however  more  complicated  than  the  single  user  case 
since  the  noise  terms  produced  are  now  not  independent. 
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The  total  cross-correlation  noise  produced  in  the  receiv¬ 


er  is  then  given  by 


K  K 


K  K 


K  K 


S2  =  E 


v  /  IV  IV  \  /  ft  XV 

53  53ei(m)en(m)}  +  E  EE  Ci(m)Cn(m)J  -2E  EE 

i=l  n=l  '  'i=l  n=l  '  H=1  n=l 

Ei(m^n(m)| 


(197) 


K  /  K  K 

=  y~](K-2i+l)PTT^m.+K2PToT-  ■-K-^T^-2-  -2E<y^  ^ 
i=l  i=l  n=l 

ei(m)^n(m)| 

The  cross  product  terms  in  equation  197  only  occur  when 

!l-m.)  fl+m  ) 

-H 


(198) 


m.  <  rn 
x  n 


The  assumed  ordering  of  the  sample  values  allows  the  cross 


product  terms  to  be  expressed  as 


T 


=  PTT 

c 


K 


E 


K 


n=i+l 


(K2-K)2PTc2 

3 


The  total  code  noise  produced  in  a  receiver  with  K  users 
can  be  found  as 


K  K  K 

-  £<K-21+1>PTTc,ni  -  £  £  PTT0(mn-mi) 

i=l  i=l  n=i+l 


+K2 3PT  T- 
c 


2  KPT  2 
c 

3 


(200) 


The  impact  of  this  code  noise  can  be  better  seen  when  the 
sample  values  are  viewed  as  random  variables.  Since  the  sam¬ 
ple  values  are  constrained  to  the  interval  [~l,+l]  ,  the 

logical  mean  for  any  probability  distribution  associated  with 
the  sample  values  is  zero.  Therefore,  the  "average"  code 
noise  produced  in  a  receiver  with  K  users  is  given  by 
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2 


(201) 


S2  =  K2PT  T-  %  KPT 
Co  c 


where  the  bar  denotes  an  average  over  the  voice  ensembles. 

The  approximation  suggests  that  the  performance  obtained 
from  multiple  codeword  modulation  formats  will  degrade  with 
large  numbers  of  users  in  conference.  The  additional  noise 
will  tend  to  degrade  performance  in  the  "large"  signal  region 
to  a  point  that  could  be  unacceptable  in  a  strong  signal  en¬ 
vironment.  The  performance  degradation  produced  is  seen  in 
the  noncoherent  bi-orthogonal  word  PWM  signal-to-noise  ratio 
approximations.  The  "large"  signal  region  approximation  with 
code  noise  from  K  users  introduced  into  the  system  is  given 
by  (Q2=PT2 /  $2 ) 

( 1-2/ft2  ) 

i.s.  (i-Kn2  +n2  (k2/2n~k/n2  ) ) 


The  "small"  signal  S/N  approximation  with  K  users  producing 
code  noise  is  given  by 


_ 1.94A2 _ 

8.54+2. 16ft2+ft2 (K2/2N-K/N2 


)  (") 


(203) 

Icoherent 


With  K<<N  the  dominant  factor  introduced  by  the  code  noise 
is  K2/2N,  the  number  of  users  squared,  divided  by  twice  the 
number  of  "chips"  per  period.  The  impact  of  code  noise  on 
system  performance  can  be  seen  by  letting  K2/2N=l/20,  which 
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corresponds  to  four  or  five  users  in  a  multiple  user  net  and 
approximately  200-250  "chips"  per  period.  The  "small"  signal 
region  expression  for  signal-to-noise  ratio  does  not  apprec¬ 
iably  change  with  the  introduction  of  the  code  noise.  The 
expression  for  the  signal-to-noise  ratio  in  a  strong  signal 
environment,  however,  is  altered  to 


a-  2/n;>  (  s\ 

(1-JsSn2  +fiz/20  \ Incoherent 


(204) 


Figure  23  demonstrates  the  impact  of  code  noise  produced  by 
multiple  users  as  the  signal  strength  is  increased.  The 
signal-to-noise  ratio  plots  show  that  a  3db  degradation  in 
system  performance  can  be  expected  when  operating  in  this 
configuration. 

The  coherent  bi-orthogonel  word  PWM  signal-to-noise 
ratio  is  also  altered  with  introduction  of  multiple  user 
code  noise.  The  coherent  system  S/N  with  K  users  pro¬ 
ducing  code  noise  is  given  by 


lN+Code  No 


ise)< 


coheren* 


_ PT2m2 _ 

2 g2  +K2PTT  -  %  KPT  2 
C  j  c 


(205) 
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Coherent 


Fig.  23.  Signal-to-Noise  Degradation  with  Multiple  User 

Code  Noise 
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In  a  strong  signal  environment  this  is  approximately  the 
same  result  seen  in  the  noncoherent  "large"  signal  expression. 
The  coherent  system  suffers  a  3db  degradation  at  £2®4.5  and 
performance  continues  to  degrade  with  increased  signal  power. 

The  exact  reduction  in  system  performance  is  a  function 
of  the  number  of  users  in  a  net.  Performance  will  not  seri¬ 
ously  degrade  with  the  introduction  of  self  noise  from  a  sig¬ 
nal  user.  However,  many  users  in  a  common  net  could  have  a 
severe  impact  on  system  performance.  All  modulation  formats 
should  be  tested  to  determine  the  exact  impact  of  this  code 
noise  phenomenon  on  pulse-width  modulation  systems. 
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VII .  Conclusions  and  Recommendations 

Conclusions 

Three  techniques  for  pulse-width  modulating  a  spread 
spectrum  carrier  have  been  explored  in  this  paper.  These 
techniques  are  forms  of  linear  modulation  and  thus  provide 
both  the  inherent  advantage  associated  with  spread  spectrum 
communications  and  the  capability  to  superimpose  several  de¬ 
sired  signals.  However,  linear  modulation  and  maximum-like¬ 
lihood  reception  produce  a  signal-to-noise  ratio  that  can  be 
improved  in  general  only  by  increasing  the  transmitted  ener¬ 
gy. 

Antipodal  pulse-width  modulation  is  the  simplest  of  the 
three  signal  formats  and  also  provides  the  best  expected  per¬ 
formance.  In  signal  space  this  format  uses  one  pseudo-noise 
code  as  a  bases  function  and  provides  a  signal  locus  of 
length  2/E  .  The  length  of  the  signal  locus  was  seen  to 

determine  the  system  performance.  Antipodal  pulse-width  mod¬ 
ulation,  however,  requires  a  phase  coherent  receiver  which 
will  significantly  increase  the  overall  system  complexity. 

Bi-orthogonal  word  pulse-width  modulation  increases  the 
signal  space  dimensions  to  provide  the  capability  to  operate 
with  a  noncoherent  receiver.  This  format  uses  two  orthogonal 
pseudonoise  codes  as  bases  function  and  provides  a  signal  lo¬ 
cus  length  of  /2e  .  All  performance  criteria  pointed  to 
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this  same  3db  degradation  in  performance  from  the  benchmark 
antipodal  format.  When  operating  in  a  strong  signal  region 
where  the  S/N>2.0,  both  the  noncoherent  and  coherent  re¬ 
ceivers  provided  approximately  the  same  performance.  The 
noncoherent  system,  however,  exhibited  a  threshold  effect,  de¬ 
fined  as  3db  below  coherent  performance,  around  S/N  =  1.5 
below  which  performance  degraded  rapidly. 

The  final  format  is  tri-orthogonal  word  pulse-width  mod¬ 
ulation  which  attempts  to  increase  the  signal  locus  length 
to  2/2E  by  adding  a  third  dimension  to  the  signal  space. 

The  effective  signal  locus  length,  however,  was  shown  to  be 
/2E  with  the  third  bases  function  only  adding  an  additional 
noise  term.  The  result  was  a  4.77  db  degradation  in  system 
performance  with  a  phase  coherent  receiver  and  severe  dis¬ 
tortion  in  the  noncoherent  version.  This  modulation  format 
has  unused  potential  that  suggests  an  alteration  could  in¬ 
crease  the  effective  signal  locus  and  allow  tri-orthogonal 
word  PWM  performance  to  surpass  the  other  formats. 

The  orthogonality  of  the  bases  functions  used  in  the 
multiple  codeword  formats  was  also  examined.  Using  the  prop¬ 
erties  of  random  binary  waveforms  the  codes  were  shown  to 
be  only  approximately  orthogonal  and  actually  introduce  a 
code  noise  into  the  receiver.  The  self  noise  introduced  by 
a  single  user  was  shown  to  have  almost  no  impact  on  system 
performance.  However,  the  effects  of  code  noise  grows  with 
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both  the  number  of  users  squared  and  the  signal  strength  at 
the  receiver.  Figure  23  demonstrates  how  as  few  as  four 
users  in  a  net  can  degrade  theoretical  system  performance. 
The  exact  impact  of  code  noise  in  a  multiple  user  net  can 
only  be  found  through  carefully  constructed  performance 
tests. 

Recommendat ions 

It  is  the  recommendation  of  the  author  that  three  as¬ 
pects  of  pulse-width  modulation  be  pursued  as  areas  of  fur¬ 
ther  study.  First,  investigate  through  simulation  or  ex¬ 
perimentation  the  impact  of  code  noise  in  multiple  user  sys¬ 
tems.  Since  theory  suggests  that  code  noise  grows  as  the 
number  of  users  squared,  a  dense  user  environment  could 
jeopardize  the  effectiveness  of  pulse-width  modulation  for¬ 
mats.  An  analysis  should  also  be  performed  to  define  the 
optimum  receiver  structure  in  the  signal  dependent  code 
noise  environment  (Froelich,  1978:1665).  Finally,  an  analy¬ 
sis  of  possible  methods  of  improving  the  tri-orthogonal  word 
PWM  format  should  be  explored.  This  analysis  should  examine 
both  methods  of  increasing  the  effect  length  of  the  signal 
locus  in  the  present  linear  format  and  the  possible  use  of  a 
twisted  modulation  formation  that  would  provide  additional 
noise  suppression  (Wozencraft,  1965:611). 
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Appendix  A 

Correlation  Properties  of  Binary  Waveforms 

The  sequences  which  are  dealt  with  in  this  paper  are  im¬ 
portant  because  of  their  correlation  properties.  The  (unnorm¬ 
alized)  correlation  between  two  sequences  of  the  same  length 
is  defined  as  (Easterling,  1964:68): 

R  =  number  of  agreements  -  number  of  disagreements  (A-l) 


The  (unnormalized)  cross-correlation  function  of  two  sequen¬ 
ces  of  the  same  length  is  the  set  of  correlation  values  of 
the  one  sequence  with  all  of  the  cyclic  permutations  of  the 
other  sequence.  The  autocorrelation  function  of  a  sequence 
is  the  set  of  correlation  values  of  the  sequence  with  all 
cyclic  permutations  of  itself. 

These  correlation  concepts  may  be  applied  to  binary 
waveforms.  In  contra-distinction  to  sequences,  binary  wave¬ 
forms  are  periodic  functions  of  time  that  extend  indefinite¬ 
ly  far  along  the  time  axis.  The  correlation  of  two  binary 
waveforms  that  have  the  same  period  T  is  defined  as: 


T 


R 


(t)C2 (t)dt 


(A-2) 


o 

The  normalized  cross-correlation  function  is  defined  as: 
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Rl2(T) 


1 

T 


X 

/ 


C 1 (t )C2 (t+T )dt 


(  A-3) 


The  normalized  autocorrelation  function  of  a  binary  wave¬ 
form  is  defined  as: 


R(t) 


1 

T 


T 


t )Ci (t+x )dt 


0 


(  A-4  ) 


The  random  binary  wave  (Papoulis,  1965:294)  which  sim¬ 
ulates  the  tossing  of  a  fair  coin  an  infinite  number  of  times 
is  the  standard  to  which  other  deterministic  codes  are  com¬ 
pared.  The  random  binary  wave  is  defined  by  the  process 
^(t)=£(t-^)  where  the  process  £(t)  is 


(1  if  heads  at  kth  toss 

(k-l)T  <t<k  T 

-1  if  tails  at  kth  toss 

(  A-5  ) 

and  the  random  variable  <g  is  uniformly  distributed  in  the 
[0,T]  interval  and  independent  of  ^(t)  .  The  autocorre¬ 

lation  function  R(t)  for  the  random  binary  wave  is  given  by 


R(t) 


(i-M),  M*t0 

0  ,  | t | >Tc 


(A-6) 


and  the  corresponding  power  spectral  density  (See  figureA-1) 
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Fig.  A-l.  Characteristics  of  Random  Binary  Wave 

(a)  Autocorrelation  Function  (b)  Spectral  Density 

is  given  by 

S(f)  =  T  sine2 ( fT  )  ( A-7) 

c  c 

The  normalized  autocorrelation  function  R(t)  for  a  pseudo¬ 
noise  binary  waveform  is  given  by 
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R(t) 


(  A-8) 


MJt|<Tc 

-  i  ,  for  the  remainder  of  the  period 

The  corresponding  power  spectral  density  for  the  pseudo¬ 
noise  wave  is  given  by 

oo 

S(f)  -  ^5-  6(f)  +  (££)£  m-)  (A-9) 

Comparing  the  power  spectral  density  characteristics 

for  a  pseudonoise  waveform  given  in  Figure  1  of  the  text 

with  the  random  binary  wave,  they  both  have  an  envelope  of 

the  same  form,  namely,  sinc2(fT)  ,  which  depends  only  on 

the  duration  T  .  The  fundamental  difference  is  that 
c 

whereas  the  random  binary  sequence  has  a  continuous  spectral 

density  characteristic,  the  corresponding  characteristic  of 

linear  maximal  sequence  consists  of  delta  functions  spaced 

1/NT  Hertz  apart, 
c 

Easterling  points  out  several  characteristics  of  the 
pseudonoise  power  spectrum  that  should  be  noted.  First,  it 
is  a  line  spectrum  with  frequencies  at  multiples  of  the 
fundamental  frequency.  Second,  because  the  binary  wave¬ 
form  is  a  constant  amplitude  square  wave  and  hence  has  con¬ 
stant  power,  there  is  a  scale  factor  inversely  proportional 
to  the  period  of  the  sequence.  Thus  if  the  period  of  the 
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sequence  is  doubled,  the  lines  in  the  spectrum  become  twice 
as  dense  and  the  power  of  each  is  reduced  by  a  factor  of  two. 
Third,  there  is  an  envelope  of  the  spectrum  which  is  deter¬ 
mined  by  the  digit  or  chip  period  of  the  waveform.  The  band¬ 
width  required  to  transmit  a  waveform  with  a  given  fidelity 
is  independent  of  the  length  of  the  sequence,  but  is  deter¬ 
mined  solely  by  the  digit  or  chip  period.  The  pseudonoise 
sequences  with  very  long  periods  are  therefore  excellent 
approximations  of  random  sequences  with  equal  chip  duration. 


The  confluent  Hypergeometric  function  iFj(a;B;±x) 
is  represented  by 


xFi(a;B;±x) 


a  (±x) 
6  1 


a  (a+1)  (±x)2 
B  (3+1)  ~TT~ 


+ 


oo 


n=0 


(a)n  (±x) 
r6)n  n 


+  . . 


(B-i; 


A  A 

where  (a)n=a(a+l) . . . (a+n-1) , (a) 0=1  .  This  function  is 

analytic  (holomorphic)  and  is  a  solution  of  Rummer 's  differ¬ 
ential  equation 


d2F 

dx2 


(6-x)  f  -o  F=0 


(B-2 


A  useful  asymptotic  development  is 


iFi(a; $;-x)= 


r(6)x 
r ( S-a) 


E 

m=0 


m !  x* 


RexX) 


r(e)x~a 

r(B-a) 


CocXoc-B-*-!) 

1 !  x  X 


(B-3 


while  functions  of  positive  and  negative  argument  are  relat 
ed  by  Rummer's  transformation 


(a; B;x)  =  eXiF! ( B-a; 8;-x)  (B-4) 

1 F i ( a ; $ ; -x )  =  e"x1F1 ( 6-a ; & ;x)  (B-5) 


For  special  values  of  a  and  6  ,  1F1  can  be  express¬ 
ed  in  terms  of  other  analytic  functions.  When  2a=B=l,2 , 3, . . . 
,iFi  is  given  in  terms  of  modified  Bessel  functions, 


1F1 (a ; 2a ; ±x)=22a-1 


T(2a) 


e±x/2 


(B-6) 
2a=l,2 .  . . 


Then  when  a  is  a  half  integer  and  B  is  a  positive  inte¬ 
ger  as  used  in  expected  value  of  the  envelope  density  func¬ 
tion,  i F i  can  be  expressed  in  terms  of  modified  Bessel 
functions  of  the  first  kind.  The  relationship  required  in 
the  expected  value  expression  is  given  by 


iFjC-JjIj-x)  =  e_x/2  [(1+x)  I0(x/2)+xI1(x/2)] 


(B-7) 
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tion  and  navigation  systems  that  are  not  possible  with  standard 
signal  formats.  The  results  of  particular  interest  are  the  abil¬ 
ity  to  screen  messages  from  eavesdroppers,  high  resolution  ranging 
and  the  rejection  of  interference. 

Three  modulation  schemes  that  exploit  these  spread  spectrum 
properties  are  analyzed  in  this  thesis.  All  three  techniques  use 
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pulse-width  modulation  to  encode  sampled  analog  voice  signals  on¬ 
to  a  digital  phase  coded  carrier.  The  performance  analysis  as¬ 
sumes  transmission  in  the  presence  of  white  Gaussian  noise  and 
CW  jamming.  The  performance  analysis  uses  maximum  likelihood  es¬ 
timation  to  define  the  sample  parameter  estimate  and  associated 
error  variance  when  probability  density  functions  can  be  defined 
for  the  receiver  sample  outputs.  When  the  density  function  is  not 
available  a  second  moment  analysis  is  used  to  define  the  expected 
signal-to-noise  ratio. 

The  results  are  extended  to  include  the  introduction  of  sig¬ 
nal  dependent  code  noise  into  the  receiver  resulting  from  the  de¬ 
gradation  of  the  orthogonal  code  word  structure.  The  impact  of 
code  noise  on  system  performance  is  seen  to  grow  with  the  number 
of  users  squared  and  the  received  signal  energy. 
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